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SYNOPSIS 


Sard's approach to optimal quadrature, resulting in a 
problem of linear best approximation, has been fruitfully 
investigated by many authors including Barnhill, ChaVla, 
Davis, Kaul, Larkin, Rabino-witz, Richter, Raina, Valentin 
and others. In the area of multi step methods for ordinary 
differential equations, however, the basic potential of 
Sard's approach, to our best knowledge, has renained un- 
explored, The present thesis may, thus, be regarded as 
a first exploratory study of Sard' s approach to optimsil 
multi step methods. 

Our adaptation of this approach to obtain 'optimal 
multi step formulae' may be summarized as follows. We 
begin with a basic multi step formula 
m ra 

1-1 J -0 

which, a priori, is known to be stable and convergent. 

We retain the a^^'s of the foimula as such and allow the 
bj * s to be variables in writing down the noam of the local 
truncation error functional at the point ^^+1* , 

the cptimal bj ' s corresponding to the point are 

obtained by minimizing this norm, (In principle one could 
have carried out the minimization also with respect to the 
a^'s. However, this over-minimization, in general. 



does not. res\alt. in a stable m\altistep method even though 
the formulae thus obtained are locally highly accurate) . 

Keeping in mind that many of the classical miiltistep 
formulae adapt appropriate quadrature formulae to the 
differential equation situation, we have also considered 
multistep formulae based on optimal quadrature. Such 
foranulae are distinct from the optimal multi step formulae 
and have been termed as 'quadrature optimal multi step 
formulae’ . 

Further, in addition to the purely optimal/ quadrature 
optimal multistep formiilae, the study also includes optimal/ 
quadrature optimal formulae subject to certain interpolatory 
constraints • 

The various formulae/methods of the thesis have been 
studied under: (a) derivation, (b) characterization, 

(c) numerical determination, (d) numerical application, 

(e) qualitative/ quantitative properties and (f) convergence 
and error analysis, 

A chapter-wise simmary of the contents of the thesis 
is given below: 

CONTENTS OF THE THESIS 

The 'body of the thesis is divided into six chapters: 
CHAPTER. 1 is devoted to' the derivation and characterization 


of quadrature optimal and optimal m\iLtistep formulae in the 
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general fraine~work of a Hilbert space possessing a repro- 
ducing kernel function. The formulae with interpolatory 
constraints have been divided into two categories: one 
with interpolatory fimctions as polynomials of a certain 
degree and the other in which these fimctions are general, 
subject to their derivatives being linearly independent on 
the node sets of the foarmulae. 

It. is shown that at each step the coefficients of the 
formxalae give rise to a deterministic system of linear 
equa'tions and that each formula is characterized by its 
being interpolatory for a certain set of functions. 

CHAPTKRS 2-5 consist of a n\xmerical impl mentation of the 

2 

ideas of Chapter 1 to the case of two Hilbert spaces H (c^) 
2 

and L (c^) of analytic functions. Quadrature optimal 
multi step formulae/interpolatory for polynomial s/general 
functions in the spaces H (c ) and L (c ) , have been 

X* JU 

considered in Chapter 2 and Chapter 4 , respectively; 
whereas, for the respective spaces, optimal multistep 
formxolae/interpolatory for polynomials/ general functions 
have been studied in Chapter 3 and Chapter 5, Each of 
these chapters contains: (1) the appropriate system of 
linear equations, particularly simplified for the space 
under consideration, for determination of multistep 

i 

coefficients, (2) the functions for which the multi step 
formulae are characteristically interpolatory (3) numerical 
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tatoles of coef ficdents of optimal multistep formulae with 
Adams ~B ash forth/ Adams -Moulton as base formulae and 
(4) tabulated numerical results when these :^l^awiLae are 
applied to a test set of twenty four differential equations, 
A numerical comparison of various multistep formulae 
occuring within each chapter has also been made. 

In each of these chapters we have studied the limiting 
behaviour of the coefficients of optimal multistep formulae^ 
as r -♦ oo • It has been established that, in each case/ 
the coefficients have limiting values equalling the 
coefficients of the corresponding usual formula which is 
of maximal polynomial precission with the same a^’s and 
other interpolatojry constraints. 

The numerical results of Chapters 2-5 show that the 
various optimal multi step formulae for the spaces H^(c^) 
and L (c^) are quite accurate and superior to the corres- 
ponding usual formxalae. The only disadvantage which the j 
optimal formulae suffer from is that -unlike in the classicaiL 
mxxltistep form-ulae -the coefficients in these do not remain/ 
invariant over the steps, / 

CHAPTER 6 deals with a theoretical error analysis of 
various multi step form-ulae derived earlier. It is shown 
that -the coefficients of the optimal multistep formxalae 
in H (c^) and L (c^) approach the coefficients of the 
corresponding usual formulae as the step size tends to zero. 
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In view of this fact, qualitative asymptotic stability 
properties of optimal multistep methods are to remain 
the same as those of the corresponding usual methods# 
However, since the coefficients of the optimal formulae 
vary from step to step, the various classical notions of 
stability and the associated stability regions, as such, 
fail to remain applicable to the optimal multistep methods. 
Consequently, for the optimal multi step methods a recourse 
is taken to the approach delineated by Henri ci. We filrst 
show that the optimal methods are convergent. Next, a bound 
has been obtained on the accumulated discretization error 
in these methods, ELnally, the propagation of round off 
errors in the application of the optimal methods has 
studied. 


been 
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quadrature optimal case 
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usual formula 
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optimal formula of a particular type 
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nxrriber of iterations in optimal type case 
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CHAPTER 1 


OPTIMAL AND QUADRATURE OPTIMAL MULTI STEP METHODS 
AND THEIR CHARACTERIZATIONS 

1,1 I nt roduct-i on 

Many of the classical m-ultistep methods are based on 
quadrature formulae with polynomial precision. Ror an 
initial value problem (I VP) 

y* = f(x,y) 

( 1 ) 

y(x^) = 

one considers 

^n^l 

yn+1” Yn-s f(x,y)dx 

^-s 

and approximates the integral occurring on the right hand 
side by a discrete quadrature formula of the type 

^ni-1 

/ f(x,y)dx = h S ^ 

n— s 

j varydng over an appropriate set. The coefficients bj's 
are choosen according to some appropriate criteria. The 
choice of the coefficients b^s in the classical quadrature 
based Adams -Bashforth/Adams-Moulton, Nystrom/Milne-Sympson 
formulae (Henrici [2?]# Jain [28], L^idus [sl] ) is so as 
to make the formula locally exact for polynomials of a 
degree as large as possible. 
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A more generalized approach is t.o consider a multistep 
method of the type 

^ yn+l-j> 

where the coefficients a.'s and b .*s are choosen so as to have 

X j 

a balance of precision and stability properties in the 
method (Henrici [27]/ Jain [28]). The precision with respect 
to the polynonial functions renders the coefficients ' s and 
b^. ‘s independent of This propearty of coefficients is 

retained even when the formulae are made interpolatory for 
the e>iponential functions e*^* However if the formulae are 
made interpolatory with respect to a general class of functions 
the coefficients a^‘s and i’j's would no more be independent 
of the point under consideration. 

Bifeginning with the work of Sard in 1949/ much work has 
been done on quadrature rules of minimum error norm corres- 
ponding to various classes of functions. Sard [42]/ Davis [ 19] # 
Valentin [44]/ Barnhill [s]# [4], [s] , [6], [s]/ [ 9 ]/ Barnhill ’ 
and Wixom [lO] / Barrar [ 13 ] and Larkin [ 32 ] have considered 
minimization of quadrature error in different manners with 
respect to weights (and nodes) to get optimal miles. Pinkus 
[ 34 ]/ Rabinowitz and Richter [35]/ [ 36 ]/ Richter [39]/ [40 ]/ 

[ 41 ] have studied propeirties of such rules. Chawla and 
Kaul [I 6 ]/ [i7]/ [is] and Kaul [ 30 ] considered such optimal 
rules subject to interpolatory conditions for ■.polynomials of 
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of degree less than or equal to the nvirriber of nodes. Kaxil [so] 
and Finney and Price Jr, [24] have studied optimal rules with 
interpolating conditions for arbitrary functions. 

No. similar approach to get optimal formulae in numerical 
solution of differential equations, however, seems to have 
been tahen up in the literature even though a need for this 
has been felt for quite some time (Gear [25] ) . 

The present thesis is an approach in this direction i 
and consists of a study of such methods which minimize the / 
quadratiore error or local truncation error in the numerical 
integration of a first order initial value probloii. 

This chapter is devoted to the derivation of the 
various formialae and their characterization. Implementa- 
tion of these formialae is done in next four chapters for 
two Hilbert: spaces of practical interest. The last chapter 
of the thesis is a study of the error analysis of these 
methods • 

We develops mioltistep methods for the initial value 
problem (1) by minimizing the norm of local truncation 
error of multi step methods in various fashicffis over a 
Hilbert space of an^alytic functions possessing a reprodu- 
cing kernel fxonction, 

A Hilbert space H of functions y(x) on a domain B 
possesses a reproducing kernel function K(x,t) if, and 
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only if, for every x 6 B the linear functional I (y) = y(x) 

iX!. 

is bounded over H* This is the well known result of 
Aronszajn (Davis [22], Theorem 12«6,1), The reproducing 
kernel function K(x/t) possesses the following property. 

For each fixed t e B, K(x,t) £ H and 

yCt) = (y(x), K(x,t)) 

for every y e H where (*,•) denotes the inner product 
in H* It then follows that the kernel function is lanigue 
and for x,t e B 

K(t,x) = KCx,t) 

In this chapter we assume that the kth derivative 
evaluation fiinctional D sf *»• f^ (x) is, for k=l, a bounded 
linear functional in H for every x ^ B. Further, we shall 
assume that H contains, as a subspace, the set of all 
polynomial functions. 

In Section 1,2 of this chapter we obtain and characterize 
multi step methods which are obtained by minimizing the quadra- 
ture error over a Hilbert space. In Section 1.3 the formxilae 
are subjected to precision for polynomials of certain degree 
alongvdth the minimization of quadrature error. In Section 
1.4 interpolatoary conditions are developed for a general set 
of preassigned fvinctions subject to which the quadrature 
error is minimized over the Hilbert space. The rest of the 
chapter follows the spirit of a general mxxLtistep method 
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where the norm of local truncation error rather than the 
local quadrature error is subjected to minimization. Thus 
in Section 1,5 we obtain the coefficients of multistep 
method by minimizing the norm of local truncation error 
over the space. In Section 1,6 the truncation error is 
minimized siabject to polyncraial precision , while in Section 
1,7 the same has been done with precision for an arbitrary 
set of preassigned functions. In all these cases the formulae 
obtained have been shown to get characterized by interpolatory 
conditions for a certain set of functions, 

1*2 Quadratxire Optimal Mxiltistep Methods (M) 

A multi step method 

lU (Si 

(2) Yn,!" Yn-s “ ' 0 < s < m, 

"" to 

where is the Kronecker delta# for the initial value 

problon (1) represents an explicit method if ts=0 and irtplicit 

method if t=l. Local truncation error T for the method 

n 

is given by 

- Yntl 

m ^ 

= - Yn_s-h S . 

^“^to 

Regarding (2) as a quadrature based multi step method the 
coefficients are ones corresponding to the quadrattare 


rule 



6 


^+1 


m 




P(x)dx^h 


iar6 


to 


The quadrature error is given by 


^n+1 m 

Q^Cp) = / PCx)dx - h E \ ^^^n+l-i^ 




iss6 


to 


and the representer of this quadrature error functional 
is given by 

(3) Q„ = ‘’I 

to 

where denotes the representer of the linear functional 

corresponding to the integral evaluation; P •* / F(x)dx ^and 

^-s 

KCx/t) is the kernel function of the Hilbert space. To 

determine the coefficients of a quadrature optimal method ( 2) 

2 

we minimize with respect to b^*s^ where 


o m _ 

(4) llQ^in = (gCt)-h S ^ 

^“\o 


Following 


m - 


g(t)-h s _J) 

j 2j 6 


to 


barkin [ 32 ] the first order c±Lange 

ro ro 

following -a change bj^, is 


6(110^11^) in 
given by 
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m - 




)) 


to 


m - 

CO 


(5) 


- h(g(t)-h S ) 

i=<5 , ^ 
to 




ro m - _ 

^="^to 


Since/ for minimizing 6(llQ^ll"^) has to vanish for an 

ro 

arbitrary choice of dOo^) , we have 


m 


(6) h S x„,.i.j,)bj = g(^»l_„) . k - 6^^(l)m 


1 * 6 , 
to 


m 


The matrix [ bhe linear sys-t&a, 

being a Gram matrix for distinct 3 * S^oCDni# is 

non~singular (Davis [ 22 ] )» Hence the system of linear 
equations given by (6) has vinique solution. 

Note that we have minimized only the quadrature error, 

l^> 

The norm of local tiruncation error T^ is given by 


m 


(7) IIT„II= llK(t.;^,^)-K(t,>,„_^)-h S I I 


where is the representer of the derivative 

functional evaluation at the point bhe purpose 

of evaluation of the norm of local truncation error it is 
assumed that derivative is a boxanded linear functional. 



Now we prove: 

Theorem 1 , The quadrature optimal multi step method (M) 
given by (2) is characterized by that it is locally Inter- 
pol atory for functions fhj(x), j = *5^Q(l)m} with 

X ^ 

hjCx) =/ K(x,Xj^^^_j)dx/ 

X 

/ standing for an indefinite integral operation. 


Proof : From (6) we observe that the multistep method (2) 
is exactly satisfied by the fxanctions whose derivatives are 

Hence the result. 


1.3 Quadrat^Ire Optimal Multistep Methods Interpolatory 
for Polvncynlals 

A quadrature multi step formiila 


C8) 


yn+l-yn-s= ‘I. 

to 


is interpolatory for polynomials upto degree p if it is 
satisfied exactly by polynomials of degree p/ that is 
by any_ linearly independent polynoiiials q^.,.<^ in the 
span fx^Jl <i <pl* Since the method is obviously exact 
for constants we donot include constants in this set. Let 


1 . (x) » TT 

PJ k=<5 




to 


^+l-k 


( 9 ) 



9 


and let 

( 10 ) 




. (x) = f l^^(x)dx # j=6^Q(l)p-l+6. 


b PJ 


■"to 


Because (x) are linearly independent and form a basis 
for span ^l/X/...,xr }/q j(x) are linearly independent 

hrJ 

and form a basis for span {x,x^, . /Xp} * Therefore we get 


m 


to 


Thus 


( 11 ) hbP = 


to 




to 


The norm of quadrature error is thus given by 

m fs 

(12) llcglU llg„(t)-h D 11 

P-^+\o' 


Let 


J”°to 




ra 


ki5p-^6. 


to 


j®=< 5 . 


to 


p-l'<'6to 


(13) s„p(t}%,^l_„)=K(t,x„^j_^)-^2:^ 


j=6. 


PJ 


to 



10 


and 

p-l+ 6^0 

- 

(14) 

■^"°to 


Then 

- 


(15) 

Ilcgu = llr„j,(t:) - 

^ to 


The fianctions ' ksp^a^^CDrn are linearly 

independent (Lemma 1«1 Kaul [ 30 I), for distinct points 
^■t*l~k' Therefore minimizing with 

respect to is equivalent to finding the best least 

square approximation to a linear combination of 

the linearly independent functions ^ ; 

k = p'**6^^(l)m}; therefore the coefficients , kssp+d^^CDm 
giving the best approximation are obtained by solving 
(Davis [ 22 ]/ Theorem (8.6,3)) the system of linear equations. 


(16) 



m 

(t)- 2 

kap+6 






1-k 


That is 


j = P + 6^Q(l)m . 


(17) 


m 
2 

ksip+ 6 


CO 

a.-.bP = 
jk k 


-"j 


to 


/ 


where we have set 


(18) ajj^= 



11 


and 


(19) 


d.. = 




Fran (18) and (19) we get 


(20) ^j]^= ®np^^+l-k'^n+l-j' "pi '"‘n+l-Jc' "np 


P“^to 

^ ^ ^ ^pi^^n+l-k^®np^^+l-i^^n+l*-j^ 

to 


and 

(21) d 


p-l+6^.o 

=^np‘=S.+l-i> - 




■j -np'”n+i-j 

*"'to 

j,k = p-*- 6^0 (Dm. 

The matrix is Gram matrix for linearly 

independent fiinctions ^ ("t?^]^) / k=p+6^Q( l)m 5 and hence 

it is non -singular. The b? /Jare •unique at all points 

Theorem 2 : The quadrature optimal mul-tistep method inter- 
polatory for polynomials of degree p<m+6^^is characterized 
by that it is locally Interpol atory for functions 

{ x^yO < i < p} U fhjjj = p+6^^(l)m J 


where 


hj(x) =/ s^p(xl3^^^_j) dx. 


Proof: Ob-viously the method (8) is exact for constants 
and {x^fl < i pJ • From (16) it follows that the method 
is exact for indefinite integral of s^^^ 
method being quadrature based. Hence the result. 
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1,4 Quadrature Optimal Mioltlst-ep Methods Interpol at ory, 
for Arbitrary Functions CM ) 

A quadrature multi step method 


m 


‘22) yn*l'Yn-s=‘' 

•^“°to 

is called interpolatory for linearly independent functions 
^ •• exact for these functions and i-f also the 
quadrature error is minimal# it is called optimal interpolator^:^ 

We assume that the matrix ^'^i ^^n+l-j i=l# j=6 

*t o 

is non— singular where denotes the derivative 

of at point • Differentiability of functions 

is a must# otherwise they cannot satisfy a 
di f f erenti al equati ons , Let 


(23) 




p-l+fito 

S * i=sl#*../P 

J-**! °to 
to 


where [c^j]= 

are linearly independent and 

p-l+<5 

^ ^pi ~ 


Obviously 


#•«,#* 


'PP 


to 


S 


C* ■ ? • J.-1 JJ 
ij J+1-& 


to 


to 


^^n+l-k^ 


X}C4-i 

The fxmctions ^^pi'*** ^ form a basis for span J 9^, , #9^1-, 

and thus if the method (22) is interpolatory for functions 
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{f 5 i-t will be interpolatory 


for fimctions 


{ 9. ,...#9 !• 

1 P 


S-ubstituting in (22) we get 


m 


All ^ -p 1L 

f (x )-f -Cx^ ) =^^i-l+6 A ^pi^^+l-j 


pi 


to 


giving 


m fo 


(25) ' 


Let denote the quadrature error in method C 22) . Then, 

proceeding similarly as in Section 1.3 we have 


iiQ^n = ngn,(t)-i^ s 






Tcs5 


to 


cu) 


in 'fs»" w p' «— 

llr^(t) - h S bj^ s (t.>Sn-l-lc 


)l» 


where 


. P-l*«to ^ 

(27) sJp(t»3^^j^_^)=*(t,x^^.l_]j5-^S, 


j-6. 


to 


to 


P-l’*^to 


(28) r^p(P) 


=g^(t)- S 


j —6 


to 


^^PJ>1-6^q' 


K(t,X^^l_j) 

P /u..“ \ , VaiD't'S. (l)m are 

Lemma 1: The functions t i^-V to 


linearly independent in H 
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Proof: Let 


m 

S 






to 


Then we define ^ linearly independent set 

of differentiable functions as in (23). Then 

1 m 


( 29 ) 


0 = (£^(t). 0) 

m _ _ 

^ to 
m 


k=p-*-6 ^ ^ np m+l-k 

^ to 

Now ^ -with [ ^i(Xn+i-j)]i^j_i non-singalar 50 IKil 




P- 


l +‘6 


to 


^mi .f. ■ ^k+l-S+.^^^n-H-k^ ^mi^^+l-j^ 

J — -t/ 


to 


to 


P-1h5, 


to 

ih~ ^ij 


6 - S 

je6 


to 


= hlc ^ 


i,k « p-f6^^(l)m . 


Substituting in ( 29) we get aj^=0, kap+6^Q(l)m, which 
proves the lertima* 

Hence the coefficients = P'**6 ^qC 1)^ which minimize 

1 1 Q^l 1 given by (26) can be determined by solving a system 
of linear equations 
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(30) 


ra 

S 




P *^F 


d!* 

J 


where we have set 


and 






Prcm these and using the properties of the kernel function 
we can easily obtain 


p- 1 + 6 , 


i -T = o 

jk np 




.•> 


to 




*pi+l-6^„^=%+l-k’ • 


to 


p-1+6 


* ®np^^+l-i'^+l~j^ 


to 




to 


j,k = • 


f* F T m 

The matrix [ajj^Jj ]^_p ^5 is gram matrix for linearly 

to 

independent functions ®nt^'^+l-k^ ' k:=p*^6^^(l)m, and hence 

^ P / 

it is non-sing\ilar* Thus bj /i can be determined ■uniquely 
from (25) and (30), Now we can prove * 

Theorott 3 : The quadrature optimal mul-tistep method (22), 
interpolating the functions ^ 9^,92# ... is characterized 
by that it is interpol a-toary for constants and 


^ * i=l(l)p } U £hj* j =sp + 
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where 

Proof : The proof is similar to that of Theorem 2. 

^ Optimal Mviltistep Methods (^) 

In this and the sxibsequent sections we develop the 
methods which minimize local truncation error instead of 
the quadrature error. Since the stability properties of a 
mill ti step method depend strongly on the choice of a^^^'s our 
major enphasis is on the development of optimal methods of 
the form 

m m A. 

(3l) ^^n+l-i"*" * 

~ to 


where are prefixed according to a known stable 

classical multistep method, the optimality criterion 
being the minimization of the norm of local truncation 
error functional. 


./v 

The local truncation error functional of an optimal 
multi step method ^plied to a function yCx) is given by 


(32) 




m m 

^n+l- , f , ^i^n^l-i- h 2 
^-1 J-o 


to 


A 

b .y* ^ 

j^n+l-j 


Therefore 
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^ rti m _ 

(33) 1 1 T^i I = S S 

±^1 J-^to 

where designates the representer of derivative 

evaluation at has been assumed that the deri- 

vative evaluation is a botanded functional in H / with 


(34) " 3x 

A 

For determining bj/j=o^Q(l)m such that the norm of the 
local truncation error is minimal in the method (3l)/we 

/V V\ 

minimize 1 mth respect to bj s* We establish: 

Lemma 2 : Let S = i KCt, 2 ^^^_^) , ^^"^'^n+l-i^ ’ ^ ~ o(l)m}^ 
for distinct points S is a 

linearly independent set. 


Proof ; We know that if i=o(l)m are m+1 distinct 

points the basic Hermite interpolation polynomials P2m'*'l/k 
Q2m^l,k; satisfying 

*ik ' ^ * 0(l)m , 


form a basis for polynomials upto degree 2(m««‘l)-l (Ralston 
and Rabinowitz [ 38 ]/ p. 70). Let 
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m __ 

(36) S + 

isO 


m 

S 

j^O 


.) s 0 


Then talcing the inner product of (36) with functions 
22rr,+l,k(^>' 

c^= 0, = 0, Ic = o(l)m, 

Completing the proof. 

Since S is linearly independent, to determine the 
optimal coefficients b. , j = 6^Q(l)m, which minimize the 
norm of local truncation error 1 IT^ I U proceeding as in 
Section 1, 2, we have 

m _ 

m 

k = 6^Q(l)m, 


Since 

] 

= It 


we have 
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m 


(37) 






to 


m 




k = 


From equations (37) we get following theorem i 

Theorem 4 : The optimal multi step method (3l) where 

a^, i = l(l)m are prefixed is characterized by that it is 

locally interpolatory for functions 

Id(x, 5„^ l-j>' J = • 

Lemma 2 indicates that if the minimization of the 
local truncation error is done also with respect to a^'s, 
the system of linear equations thus obtained will have a 
non~singular matrix and that such a formula will also be 
unique. However such formtilae may not possess the stability 
properties required in a multi step method for the whole 
internal. Such a formula may be useful where the major 
interest is the minimization of the local tfuncation error 
( e, g. in the beginnir^ of interval or at the point where the 
solution may be fast changing) t nevertheless. 

If the minimization in (33) is done with respect to both 
a. * s and b.‘s # the normal equations are 

^ J 
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m 


m /V 


1=6 


to 


, K(t^ ^ i=l(l)m 


and 


m 


E (D(t,x^ ^_j), K(t,x„ j_^))a^ 
k=l -' 


m 




to 


(D(t,5^^1_j), K(t,5^^^)) , j 


= ®to‘l>"' 


which simplify to 

m 
2 

kal 

(38) 


I" - ,/v 


Is 6 


to 




m 


m 


l=«to 


J " ®to 


( l)m. 


From these 2m equations we conclude: 

Theorem 5 ; An optimal multi step method given by 


m 




m A 


‘39) yn*l= Vn*l-i* 1-J' yn,i_j)- 


to 
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in which the minimization is done both with respect to a^'s 
and is completely characterized by that it is 

locally interpolatory for the set of functions 

; i=l(l)m} U ; j ¥= 6^^(l)m } 

1,6 Optimal Multi step Methods Intejrpolatorv for 
Polynomials (M^) 

For the purpose of derivation, a multi step method 
m m 

yn+l = f 

*“ ^ to 

can be written^ using divided difference notation# as 

m m 

^ to 

If p ^ m we can fix { j = d^^CDp-l + “^to^ ^ manner 

that the method is exact for polynonial upto degree p. 

Thus we have a formula 

m P"’^'*’'^to 

x-x J"°to 

m 

J=P»«to 

in which the minimization of the norm of error is done 
with respect to J = 

f 

If D(t^x) is the representer of the derivative evaluation, 
which is asstmed to exist, then the representer of the 
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divided difference Y* obviously 

to 

to “* 

which is obtained by having the divided difference of D(t/x) 
with respect to the second variable. This being simply a 
linear canbination of ^ involving 

constants# for any function f(t) we have 


(f (t) , D [tjx 


n+l-d 


|r • # 


to 


h+l 


-P> 




and therefore 




# • • • # 


to 


%i+l J ' ^ 


(42) 


= D*[ 


X, 


•n+l-6^^'"’'*'^'^l-j 

to -* 




'•n-^l-i 


-1^ ' 


which is the Ci/j)th entry of the matrix obtained in the 
minimization of the norm of error in method given by (41) • 
Thus we can write the normal equations as 

(43) Px = q 

where P is an (m-p-d^^+l) x (m-p-6^^+1) matrix where 
p. . are given by (42) and 

^ J 


% = 




■^n+l' • * • '^+l-i 


( 44 ) 


m 


"kf 1 * ’ * ^^+l-i^ ^ 

S ^n^l -fill*" ^•*'1-1^ 


Istfi . 


'to 


to 
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The matrix P is obtained by performing elementary row and 
colxmn operations on a leading submatrix of the matrix of 
the linear system (37). Therefore P must be non-singral ar 
and thus the method (41) is also uniquely determined. By 
writing the representer of local truncation error for 
method (41) we can easily prove J 

Theorem 6 : A multistep method given by 


m ni /V. 

to 


is optimal interpolatory for polynomials of degree p < m 
if and only if it is interpolatory for the functions 
{x^jk=0(l)p} U (hj# j=p+<5^^(l)m} / where 


= d[x;S^ j_S . 

to 


x 


•n-6. 
to 


J 


1.7 Optimal Multi step Methods Interpolatory for Arbitrary 
Functions (M ) 

In this section we give a different approach for deter- 
mining and characterizing a mxiltistep method interpolatory 
for a set of arbitrary functions, (As compared to Section 1.4) . 
The approach is based on Lagrangian mxaltipliers. We £Lrst 
prove a result of general applicability: 

Iheorom 7 : A necessary and sufficient condition that the 
multistep method given by 

m m 

Vn+l-l'' *' 


I 
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be interpol atory for functions is that for 

every 9* e span { f^^ , . f^ } for which 

there holds 

m 

’’ '■^*1 > ’’‘^♦ 1 - 1 ^ 

Proof : The necessity paxt is obvious. For sufficiency 
assume that the stated conditions hold. Consider the system 


(45) 


m 

S 

jiSsS 


m 




to 


of linear equations in bj * Applying elonentary row operations 
on the system (45) we can reduce it to one in which the 
matrix is in a row reduced echelon form. In this system the 
existence of a non -zero augmented entry corresponding to a 
zero row of the matrix would contradict the stated conditions. 
It follows that the system is consistent , conpleting the 
proof of the theorem , 


In view of the above theorem and its proof# if a multi- 
step method interpolatory for span {f^/.,,,fp} is possible# 
there_exist q (m+d^^^) functions g^/g 2 '***' 5 q 
span If^#,..#f } such that 


rs«l^ Ogi^x^+l-.j)]i=l(l)q ^cFdira[span 

J=6^^(l)m 


( 46 ) 
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and that the multistep method is interpol atory for 
span } if^ and only if, it is interpolatory for 

a. 

Hence with preprocessing, if necessary, we 
can assume that fj^,,,.,fp themselves are such that 

(47) 

Here in case p = m all the coefficients j ‘ s get 

uniquely determined and no optimization of error will be 
possible. Thus for optimal interpolatory methods we will 
assime that p < alongwith (47). 


For an optimal method interpolatory for f^ span {f-,.,,,f } 

Jr 1 ' p 

we have 

llT^ I 1 is minimum 
n 

subject to 

(48) (T^, fj^) = 0 , k = l(l)p, 

A.p 

where T_ is the local truncation error functional corres- 
n 

ponding to the method. 

Equivalently we minimize 


(49) (?^. $^) ♦ 

where ^'s are complex Lagrangian multipliers. Thus 
we have to minimize the (quantity 
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m 


m 




to 


x-x -^“^to 


m 


♦ ^SACK(t,x„^l)- s 
k=l i=l 


in - 


J-o. 


to 


giving rise to a system of linear equations 


'A # 


'b" 


~c “ 

F 




f 


(50) 


where A is an (”''*•6^2^) x (m + matrix with (i,j))th 

entry ^ is p x (m-«-6^^) matrix with 

(k,j)the entry as adjoint of F, 

0 is a p X p null matrix , 

b = [b|^„....,b^]^ X =[\l f = [fi ^ 


m 


where 


x^l 


and c = Cca^o#- • * where 


m 


a = i-j> • 


If the rows of P are assumed to be linearly independent 
from that of A the system of linear equations can be solved 
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to have a unique solution, thus obtained vd.ll satisfy 

the required optimal interpolatory conditions. Now we 
characterize the optimal interpolatory multi step method M . 


For simplicity in the form of characterization we will 
assume that the last p colvamns in the matrix in (47) are 
linearly independent. (Note that in the general case sane 
p coltomns of the matrix are linearly independent and in the 
corresponding characterization the indies as of these columns 
instead of the last p columns in the simpler case would 
appear in the formiilation.) 

Let the interpolatory conditions appearing in (50) be 
rewritten as 

[eIp] [gl = f 


where 






J=sim+l-p( Dm 




F 

'm-p+1' • 


. f is 


as in (50). We get 


e * P’'^(f - Ed) . 


Let 


(51) 


r a P^^f 



28 


and 

( 52 ) 

We have 

( 53 ) 


G = P 


e = r - Gd . 


Then we can write 


m m-p ^ , 


■■ ^n+1 


m 

m-p 

Vfn*X-L- 

2 

J^to 

m 

m— p 

jssm-p^l 1^6^^ 

m 

m 


1? 


'to 




m-p yv 


^ } 


jr p ^ 1 *1 

" ’"j Is ’"J ■^'^' C 

j to ^ 


With the notations and assumptions as in the above the 
optimal interpol atory multistep method is characterized as 
follows i 

Theorem 8: A multi step method given by 


^n+l “ .f. Vntl-i*^ ^.^5 '^n+l-j^ 

i-x -* to 

where the e^‘e are prefixed^ is optimal interpol atory for 



29 


the functions of (47) if/ and only if# it is 

interpolatory for the functions 


Cf^,.../fp} U £h^. ; j=: 5^^(l)m-p} # 


where for j = 6^^ ( 1 ) m~p , 


m 




-P^l 


D(x,z) is the representer of the derivative evaluation 
fianctional at point z and g's are as in (52), 


Proof ; Equation (54). shows that to determine optimal 
coefficients (bj/ j=5pQ(l)m— p) / the normal equations 
may be obtained by minimizing the representer of (54) with 
respect to the coefficients bj/ j=6pQ(l)m-p, 

Thus we obtain the equations 

(hj, T^) =0 , J = 6^.o(l)m-p, 

where 

m 

hj - ♦ 2 5i.^4p, • 

ism-p-cl 

As the method has been made interpolatoary for functions 
I f^,f 2 #«**/fpJ the proof is complete* 

In Chapters 2 to 5 we have implemented the methods 
of this chapter in different situations* We present there 
the coefficients for 5 step methods with s = 0 for the 
quadratiire optimal case and a 2 = a^= a^s 0 for 
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optimal case. Coefficients for non~interpol atory situation, 
interpolatory for polynomials upto degree 1 to 5 and inter- 
polatory for e:5^Cax), a = ± 1.6 have been calculated and 
presented in the form of tables for quadrature optimal and 
optimal case in the respective Hilbert spaces. Coefficients 
have been obtained for the interval [-2.0, 2.o] with step 
size h = .1 and equispaced points. All the methods have 
been applied for a set of twenty four differential equa-cions 
by providing exact starting values (calculated frcm the 
analytical solutions) at 5 points prior to -1.2 . The 
solutions have been computed on the interval [-1,2, 1*2] 
and ccffnpared with usual 5 step methods. 

The twenty four equations and the analytical (particular) 
solutions are as follows: 

1. g = xy3 ; y = (2.0 i2- 

2. g = xy3 J y = ( 2 . ll^-x ^) “^2 

3. ^ 2xy^; y = ( 2.03^-x^) 

4. ^ = 2xy^; y « (2. 11^-x^) 

5. ^ = -2x.e~^f y = log( 2.01^-x^) 

6. ^ = 4xV^ t Y = x^)"”^ 

g = 6 x\2 , y ^ (2.01^- x")“^ 


7. 
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• y = (2.01'^- 

3x^Y^ J Y = (2.01^- x^) 


= 2 xy * y = 


= 3(y-“^^) # y = cosh (3x) 

= sCy^-e^^) ; y = sinh (3x) 

= XY^ > y = (2.75^-x^)“^/^ 

= xy^‘ ; y = (2.85 ^~x2)“^/^ 
= 2xy^ j y = (2.77^- x^) 

= 23?y^ ; y = (2,85^~x^) ^ 
a 6x^y^ J y = (2,85^-x^) 


3x y 


5 3 . y ^ C2.9 "-x")-^/2 


^ 2x(y2- ^ ); y = (2.022 -x2)-^+ (x 2 ^. 9 )"^ 


2 „ 

X ^*.9 


2xy^( 2x^-3. 1804) ; y = (2.02^-x^)"^ (x^+.9) 
xy^(2x^-3.1804) 5 y = (2.02^-x^) Cx+.9) 

X C(y-' (' • ^ ■■ " ■) y=( 2 . 02 ^-x^)"^^^+(x^+. 9 ) 

'^x^+-9 '^x^+,9 


x^+.9 


2 rv 

X •!■ ,9 


2 2 
,2*02 -Xx 

; y a log ( 5 ) 


24 . ^ 

dx 


ci Q . 2 .2 2xv C 2.0 2 -X 

3x y-^(x^-*-.9) - t Y - ^—x 


6 6. -3-/2 


2: ^ 

X +.9 


x^-«-.9 
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The first nine of the solutions have singularities 
of various types near the boundary of the domain of the 
Hilbert spaces that would be considered in chapters 2-5, 

The solutions of next -tSaree equations 10-12 are entire 
fixnctions of rather fast growth. The solutions corresponding 
to equations 13-18 have singularities comparatively far 
remote from the domain of the underlying spaces. The last 
six equations have solutions with singularities inside the 
domain. Thus through these twenty four cases we have tried 
to cover most of the important theoretical situations that 
arise in the case of analytic solutions. 



CHAPTER 2 


QUADRATURE OPTIMAL MULTI STEP METHODS IN H^(c^) 

2, 1 I nt roducti on 

Quadrature formulae of various types for functions in 
the space H'^Cc^) have been studied by Larkin [32], 

Richter [39], [40], [4l], Chawla and Kaul [l6], [l7], [is], 

Finney and Price Jr [24] and others. The methods involve 
choosing the coefficients and/or nodes to minimize the 
noirni of error without or with some constraints* The resrilts 
of these authors show that the error in the quadrature 
form\ilae optimized over H (c^) is substantially smaller 
than that corresponding to classical quadrature rules. 

In this chapter we consider the applicability of such 
optimal quadrature rules in H (c^) to initial value 
problen for first order differential equations* 

For the purpose of numerical illustration we have 
worked with quadrature optimal multi step methods of Adams - 
Bash forth and Adams -Mota ton type. 

2 

In Section 2.2 we describe the space H (c^) and 
establish the boundedness of derivative evaluation 
fxanctional of any order at a point inside c^* In Section 
2.3 we present the construction of quadrature optimal 
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2 

methods in H (c^) and some nxjmerical results obtained 
on applying the method of S-step on differential equations. 
The coefficients and norms of local truncation error 
functionals of implicit and e^^licit 5 step methods for 
h = 0.1 and the results obtained by applying these to 
differential equations are given and compared with the 
corresponding usual metnods. In Section 2,4 we describe 
the quadrature optimal multi step methO(^ of 5 steps inter- 
polatory for polynomials of degree one to five. The 
coefficients and local truncation error functional norms 
are given in tabular form and behaviour on the differential 
equations is described. In Section 2.5 we present quadrattire 
optimal multi step methods of S—step interpolatory for a 
set of selected functions. In Section 2.6 we discuss the 
behaviour of the coefficients of quadrature optimal miolti— 
step methods in H (c^) as r -* <» • In Section 2, 7 a 
con^arison of different methods is presented, 

2,2 The Hilbert Space H^(c^) 

2 

The space H (c^) consists of analytic functions 
regular in the disk = {z i lz| < r} with the magnitude 
of boundary values square integrable on the circle 
Cr = {z: Izl = r}. The space H (c^,) is a Hilbert space 
with the inner product defined by 

(f#g) -f f(z) g(z) ds, f,g e H^Cc^, 


( 1 ) 
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the integration being with respect to the length element 
on the circle. The space possesses a reproducing kernel 
function given by 


(2) KCz^t) = ^ 

2t(r - zt) 

and a complete orthonormal sequence of functions ([43]) 

(3) ^vCz) = - - r- . k=0,l,2,... 

^ v/‘2ur r^ 

2 

By virtue of the space H (c^) possessing a repro- • 
ducing kernel function the point evaluation functional 

L : f(eH^(c )) •*> f (z), (z e D ) / 

Z XT IT 

2 

is a bounded linear functional on H (c^) , In the sequel 
we require the boundedness of the derivative evaluation 
functional of order k 

: f(eH^(c„)) f (z) (z e D ) , 

^ X XT 

where k is a positive integer: 

Theorem 1 : For z e and k any positive integer 

k 2 

D is a bounded linear functional in H (c^) with 

Z XT 


■ C5) 


I ID^(z) II 


•« '■•tar ICai* ^ ^ 

'/27V.(r- 1 z I ) 


Proof : Using Cauchy integral formula# if 6 denotes the 
distance of a point z e from c^# we have 
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I £« (z)I < i / 

r 


Idw I 


k! 


271 0 c. 


!f(w) 1 Idw I 


< (£— ) ^ ^ • 1 1 f I I 

- '•271;^ ^k+1 

O 

by Cauchy-Schwarz ineq;u.ality. From this the result is 
immedi ate . 

That the operation of evaluating a definite integers! 

2 

corresponds to a bounded linear fxanctional in H (c ) is 


well known (e.g, Larkin [ 32 J) moreover we have 

( 6 ) 


^2 r ^ — tz 

/ f(z)dz = (f(t)/ ^))#feH^(cp 

Zi r - tZ 2 


showing that the representer 1/ _ v (t) of the definite 

ji Z 2^ 

integral fvinctional is given by 

r^— tz 


2.3 Quadrature Optimal Mxiltistep Method in H^Cc^ 

For determining the coefficients of a guadrature 
optimal multi step formula 


( 8 ) 




2 

in H (c^) we , use the equations (1.6) (equation (6) of 
Chapter 1)* (2) and (7), Thus we have the normal equations 
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(9) 


ro CO <s> 

C b = d 


with 




and 




27ih X 




ro 

d. 

1 


log Vs , 3 . tkl ),lf ^ 0 


=S +1 


X 




2%rh. 




to 


Equations given by (9) can be solved analytically as follows: 
We know that 


(10) det (C) = r 


where 




m 

TT 


m 






and 


k = (ra+6^^ ) (m+6^^-1) ^ 


M = (|^) 


tl 


Using Cramei^s riile we get c”^ = g# (say)# where for l,k=6^Q(l)m 
m „ m o 


\ “I 

%k“ ^ 2 Tl 


1-A j 1-5 


1 = 6 . 


to 


■aTrn^d^ 2 


to 


m 




J=A 


j'A 

m 

T 

i=6 


to 


TT ^^^l-l”^n+l-i^ ^ 


i/^k 


to 
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There-fore. the coefficients are given by 


- 1 r f 

h “ h 3lk 2nh * 


-L n / n^l-kx ^ ^ 

log (-3 )i / l=6^^Cl)m, 

^ntl-k ^ “^ntAtl-k 


{ o denoting limiting value if 


0 


Thus the method (8) gets determined e 3 <plicitly. 

In view of Theorem 1.1 (Theorem p q denotes qth theorem 
of fSiapter p) we have 


Theorem 2: The quadrature optimal multistep formula (8) 
2 

in H (c^) is characterized by that it is locally Inter- 
pol atory for fxinctions 


r 


yj^(x) 


log(r -X . if ^ 0 


X , if x^^^„^=0 . i=6^^(l)m. 


Proof ; Fran ( 2) » (7) and Theorem 1,1 it follows that the 
formxila (8) is interpol atory for functions 

r 


g^(x) 


f ^ ^n-i*l-i 


log(r if 0 


X 

27cr 


* ^Si^l— i “ ® ^ • 


Suppressing the scalar mialtipliers in g. (x) we get the 

wLl 

result. 

From this theorem we note that if 0 for 

some i with 6^^ < i < nv, i.e. the origin is a node used 
in the quadrature optimal multi step formula, it becanes 
exact for polynomials of degree 1. 
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In Table 2,1 we present the norms of local truncation 
error functionals and the coefficients of esq^licit and 
implicit quadrature optimal 5-step method in H (c^) with 
s=0. The norm of local truncation error for 5-step usual 
methods is also given at each point. The norms of local 
trruncation error functionals as defined in equation (1.7) 
(equation 7 of Chapter 1) may be calculated by using the 
following: 


For any coefficients b^, j=6^Q(l)m the norm of the 
local truncation error functional T^ at is given by 


2 r 


"^n'‘ = 2t 


[ 




- 2 Re(- 


*) 


r 


'^■hl =^n-s 


(11) 


m 


m 


- h S bj 2 


where 


m 


- h s b, - J] , 


1 . 4 O 


^n-^l-j^n+l-h^ 


g(x) = x[(r'^-5^^^x) 


-2 


Here Re( ' ) denotes the real part of the quantity. 

All the computations were carried out on DEC-1090 
computer system of the Institute. The coefficients (which 
have been obtained by Gaussian method for the system of 
linear equations) and the nornns of the error functionals 



40 


were computed using double precision, while the solutions 
of the differential equations were obtained using single 
precision. All through the thesis the computations are 
confined to the interval [-2,2], a real situation. We 
have tahen r = 2.01^h = 0,1. The interval on which any 
formula is applied is [-1.2, 1.2](by providing true values 
of the solution at the five points —1.7, -1.6, —1.5, -1,4 
and -1.3)» Implicit formulae in all cases have been 
implemented by providing initial approximation generated 
by an application of the explicit fomula of corresponding 
type on the last five values obtained by the implicit 
method (or the exact values at the beginning). Iterations 
were carried out till the difference of two consecutive 
iterates was smaller than lx 10 in absolute value, 
s\xbject to a maximxm of ten iterations. 

From Table 2.1 we observe, in both implicit and 
explicit cases, that at x = —1,5 the norm of local 
truncation error functional for usual formiila is more than 
a thousand times the corresponding norm for quadrature 
optimal formula. In central region (around x = 0.2, 0,3) 
the norm of local truncation error functionals of the 
usual and quadrature optimal formulae are almost same. 

Prom here onwards again the error functional norms for 
the usual formulae start beccmming larger than those of 
the quadrature optimal formulae. 
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From Table 2»1 we also see that the coefficients of 

2 

the H ( c^) —quadrature optimal methods vary from point to 
point. It is interesting to note that the magnitude of 
each of the coefficients/ however, is an increasing fxanction 
of X except the first coefficient of implicit method which 
decreases with x. Moreover in the central region the 
coefficients approximately equal the usual coefficients. 

Following Table 2.1, in Tables 2, 1(a) --2. 1(d) we 
present numerical resxxLts obtained by using the coeffici- 
ents of Table 2,1 and the usual coefficients on the twenty 
four differential equations (given at the end of Chapter 1) , 

2 

From Tables 2, 1(a) -2, 1(d) it is clear that the H (c^)- 
quadrature optimal 5-step methods have performed better 
than the corresponding usual methods which txxm out to be 
the Adams-Bashforth/Adams-Moulton 5-step methods (for s=0) ♦ 

I The explicit method has been just better on equations 
- 13/ 14 and 24 one decimal place more accurate on equations 
2# 3# 6/8/9/ 10/ 11/ 15/ 16/ 17/ 18/ 19/ and 23 ; and two decimal 
places on equations 1/4/5/7 ,12/20,21 and 22, The implicit 
quadrature optimal method is foimd Just better on equations 
17 and 18; caae decimal place better on equations 7,11/12/13/ 
14,15/16/19/20/21,22/23 and 24> and two decimal places on 
equations 1, 2/3/4/5/6/8/9/ and 10, We also observe that 
the full 10 iterations were carried out rarely. Indeed, in 
this work idea has been only to check the effectiveness 
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of the basic implicit scheme. In actual practice^ of 
course# one wovild utilize a more appropriate predictor 
corrector scheme. 


2 

2.4 Quadrature Optimal Multi step Methods in H (c^) 
Interpol atory for Polynomials 

In Chapter 1 we have giyen an approach for deter- 
mining quadrature optimal miiltistep methods interpolatory 
for polynomials in a general set up. Here we specialize 
the case when the points along a straight 

line in the consecutive points eguispaced with distance 

h. It is obvious that 

(12) y^^i = h yn^6^i^ 


represents a general m step method of quadratixre type* 
Frcm this it follows that a general method interpolatory 
for polynomials of degree p < m +■ may be written as 


(13) '’s ^ / Y 


jsO 

m-6, 




to 


= h S f( 


X. 


J=P 


-■tl 


■n+6. T '^n+'d. 


tl 


where the coefficients ')j* s are defined as in the 
corresponding usual formixla. In a (quadrature 
optimal miiltistep formula interpolatory for polynomials 
of degree p# the remaining coefficients j=p(l)m-6^Q 
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are to be determined by minimization of the norm of 
quadrature error functional over the space H KC^} . 
Proceeding as in Section 1.6 it follows that the unknown 
r ' s are to satisfy the normal equations 

j 


(14) 

where 


oP P = 'S’ , 


= ‘’'P ' 

^3 2 7 C 1=^0 


.“I 




^n4.6^,-l 




^ ( 1 s (-i)^’*’^(i)(K^ 


^ k:s0 1=0 




i/j = pCDm-S^o- 


Here replaces i. J* If " ° ' 

The .multi st«> ooefficlehts b. (as in (1.8)) can then be 

calculated as in the \asual case. 

Theorem 3! In H^(Cj.),a quadrature optimal multistep method 
(13) Interpolatory for polynomials of degree p c m+6^.j^ rs 
dharacterlsed by that it is locally interpolatory for 
functions (A 1< j <P>U (hjW.- J = pCD^-^to^ 
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h.Cx) = i log (r^-x J J }* , 

' =S,.6^,-k 

where * indicates that if 0/ the expression 

within the brackets is to be replaced by x/r^. 

Proof ; The result is immediate from (14) along the lines 
of Theorem 2» 

In Tables 2.2 through 2.6 we present the coefficients 

2 

and error norms for H (c^) -quadrature optimal 5— step 
methods interpolatoiry for polynomials of degree p = 1(1)5, 
For the comparison purpose# error norms of the corresponding 
usual methods have also been included. While going through 
Tables 2.2 to 2. 6 we observe that the error norm for qua- 
drature optimal method starts increasing when the interpo- 
latory conditions are increased. In Table 2,6 # wheret the 
polyncKTiial precision is 5, the explicit method has become 
the usual method. 

As the quadrature optimal m\altistep methods inter- 
polatory for polynonials of decree p = 1(1)5 have performed 
more or less in between the optimal and the usual methods 
we do not include complete tables of the results obtained 
by applying these methods on the differential equations. 
However a simmary of the restxlts observed is as follows ; 

The H^(cp -explicit quadrature optimal multi step 
method with p»l has been just better on equations 13# 14, 
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16 and 241 decimal place better on equations 2# 3# 6,8/ 

9/10,11,12,15,17,18/19 and. 2 . 3 } and two places better on 
equations 1,4, 5, 7, 20, 21,aind 22.Th.e implicit method -with 
p=l is found just better on equations 17 and 18; one decimal 
place better on equations 7,. 11, 12, I 3 , 14, 15, 16, 19, 20, 21, 22, 23 
and 24; and two decimal places better on equations 1,2, 3,4, 
5, 6,8,9 and 10. 

The explicit method with p=2 has been just better on 
equation 24} one decimal place better on equations 2, 3,6,8, 
9,10,11,12,13,14,16,17,18,19 and 23} and two decimal places 
better on equations 1,4,5,7,15,20,21 and 22* The implicit 
method with p=2 was just better on equations 17 artd 18> 
one decimal place better on equations 1,5,7,11,12,13,14,19, 
20 , 21 , 22,23 and 24} and two decimal places better on 
equations 2,3,4,6,8,9,10,15 and 16, 

The explicit method with p=3 is foxand just better on 
equations 11,12,14,16,17,18 and 24} one decimal place better 
on equations 1,2,3,4,5,6,7,8,9,10,13,15,19,20 and 23} and 
two decimal places better on equation 21. The implicit 
method with p= 3 has been just better on equations 11,12, 
17,18 and 24} one decimal place better on equations 1,3,5, 
7,13,14,15,16,19 , 20 , 21,22 and 23} and two decimal places 
better on equations 2, 4, 6, 8, 9 and 10. 

The explicit method with p=4 was just better on 
equations 6 , 14 , 16 , 17/18, and 24} one dedmal .place better on 
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ecjuations 1,2,3,5,7,8,9,10^13,15,19,20,21, and 22j and 
two decimal places better on equations 4 and 23. The 
implicit method with p=4 has been just better on equations H, 
12,17 and 18?one decimal place better on equations 1,2, 3,4, 
5,6,7,8,10,13,14,15,16,19,20,21,22,23 and 24? and two 
decimal places better on equation 9, 

The implicit method with p=5 has been just better on 
equations 3,10,13,15,17 and 18; one decimal place better on 
equations 1,4,5,6,7,8,9,14,16,19, 20,21,22,23 and 24; and 
two decimal places better on equation 2. 

The e3{plicit method with p=4, however, has gone just 
worse for equations 11 and 12. The same has been observed 
in the case of in^licit method with p=5, 

2 

2,5 Quadrature Optimal Multi step Methods in H (c^) 

Interpol atory for a set of Preassigned Functions 

The coefficients bj of a quadrature optimal multi— 

2 

step formula, in H (Uj.) , 

yntl- yn-s= I’ yn+l-j)' 

interpolatory for linearly independent functions ...,9^} 

can be obtained by solving the system of linear equations 


" 5 

F 


-gi- 


3 

F 

1m» 

0 


_x_ 


e 
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where C and d are as in equation (9), / 


^to 




«to' 




and 0 is p x p null matrix* 

Following theorem gives characterization of a 

p 

quadrature optimal mxiltistep method in H'Cc^) inteirpola- 
tory for a set of preassigned functions. 

2 

Theorem 4; A quadrature optimal multi step method in H-Cb^ 
interpolatory for functions {9^^, , , , ,9^1 is characterized 
by that it is locally interpolatory for constants and 
functions 

U {h^(x), k=p+6^Q(l)m } 

where 

hj^(x) log 

p-'l+S^o 
S { S 

j^6to q*«5to 

. f log .x) }* 

where 

w =* Cwijli^KDp 

j=6to^l)P~l‘^\o 


'^j+‘l~5pQ/q *^q+l-6^^^^+l-lc^ ^ 
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wi t±L 

^ “^to ^ ^ ^ 

and {. }* means the expression to be replaced by ^ if the 

r*^ 

denominator in the expression is zero. 

Proof i The proof follows from Theorem 1,3 after a little 
simpli fi cati on. 

Table 2,7 corresponds to the quadrature optimal case 
interpol atory for functions ‘P^(x) = e>p(l,6x) and “ 

exp(-^l,6x). The behaviour of the error norms and the 
coefficients is fo-und to be as in the case of quadrature 
optimal methods. The explicit method was found just better 
on equations 17# 18 and 24; one decimal place better on 
equations 1, 2#3#5#6,7#8#9# 10#11#12#13#14#19, 20 and 23? and 
two decimal places better on equations 4#15#16#21 and 22 as 
compared with the usual method. The implicit method has 
been just bett<=r on equations 17# 18 and 24; one decimal 
place better cxi equations 1#4#5#6/7#11#12#13#14#15#16,19#11# 
20# 21# 23# and two decimal places better on equations 2# 3# 
8#9# and 10, 

2, 6 Limiting Behaviour of the Coefficients as r 

Frcxn Tables 2,1~2,6 we observe that the coefficients 
of the quadrature optimal mtiltistep formiola/interpol atory 
for polynomials of a given degree are approximately the 
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same as the corresponding coefficients in the usual formula 
for the points under consideration near the centre of the 
circle c^ (or loosely speaking when the points are well 
inside the domain of the space) . Properties of coefficients 
such as this have been considered in [?]# [24]/ [ 30 ], [39 ] 
and [40]. For the space H^Cc^ / Kaul [so] and Finney and 
Price Jr, [24] have studied the behaviour of coefficients 
in quadrature formulae as r , The basic approach for 

such a study is given by Barnhill [?] , 

The following results give analogoues of such beha- 
viour of coefficients for quadrature multistep formulae 
in H (c^) interpolatory for polynomials/ arbitrary functions 
as r •* «> , The basic tool in the analysis is provided by: 

Theorem 5: Let 

^ F 

(16) Yntl = .f, ^ yn+l-j> 

J”“tO 

be a quadrature type multistep formula interpolatory for 
and let the matrix L j- 5 ^ 

be non- singular-. Let 

I 

^n-cl ^n-s"*" ^ ^jr ^^^n+l-j'^n+l-j^ 

2 

be act . quadrature optimal multi step formula in H (c^ 
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interpolatory for >/ 0 <P < 

(p=0 implies non -interpolatory case for these f\mctions). 
If the quadrature error functional for (17) satisfies 


(13) 

lim ^(9!.) = 0 , J = p+6 (Dm , 

then 

lim = b^ / j = 6, (Dm* 

x-* 00 -* 

Proof: 

Let 






to 

and 

0,= W L 


Since (16) is exact for the functions 9^ ^...#9 and 

to ™ 

(17) is exacl for 9^ /.../9^ - . * , writing the 

formulae (16) and (17) for the functions 
and subtracting we get 

hGB^ + hGB^ = 0, 

In view of (18) and the interpolatory nature of (17) for 

the functions ^ ~ ^ to^ ^ 

r -*> 00 , Hence 
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lim hG(B^ - B^) =0 
r-* oo 

Since the matrix G is non -singular we have 

lim B^ = B , 

r ' 

r-* oo 

which proves the result. 

Theorem 6. In H^(c^) the coefficients of a quadra- 

ture optimal multi step formula 

m 

(19) f h s 

^“®to 

approach the coefficients of the corresponding usual 
formula as r -* oo . 

Proof : Let denote the quadrature error functional 
of a usual formula 


( 20 ) 


m 

^n-s"*" ^nfl-j^ 

3=6 to 


Tlhen 

=0 , 0 < i^ • 

In view of Theorem 5 we need to prove that the quadrature 
error in (19) satisfies 

lim \(x^) ==0# 0 < i • 

r 

since (19) is quadrature .optimal/ I 1 0^^ I 1 < * Hence as 

i n as defined in (3) is a complete orthormal 

Km iC^sQ 

secguence in H^(c^)/we have# equivalently. 
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i=^0 ^ ^ k=0 ^ ^ 


giving 


“ -2i , is ,2 ~ -2k / k. .2 

S r ^iQ (x^)l^ < S r IQ (x ) 1 . 

i:stO k=m+l-6^^ 

to 


Taking a fixed i# 0 < i ^ m - 6^^ we get 


lQ (x^)!^< S IQ (x^) 1^ 

- k^m^l-6^^ 


in which 


the expression on the right side approaches zero as 
r -»• o® t completing the proof. 

Theoram 7: In H^(c^)# the coefficients bj^ of a 
quadrature optimal multi step formula 


( 21 ) 


^n+l' 


!= y + h 
^n-s 


ra 

S 

j-6 


jr 
to 




■n+l-j 


) 


interpolatory for polynonials of degree p < ni‘*‘l“^to * 
approach the coefficients of the corresponding usual 
formiola as r *♦ <» . 

Proofs The quadrature error fxinctional for ( 21) 

satisfies 


(22) =0 / 0 < i < p-1 

<N>P 

Proceeding as in the proof of Theorem 6^ < 110^^11 


gxves 


S r 

iap 


~ k=inipl- 


•2k 


Ic X 
lQj,(x^)l. 


to 



For any value of p < i < we have ^ therefore/ 

k=m*l-6^„ 


Thus 
( 23) 


. / i ^ 

lim Q„(x ) = 0/ p < 1 < . 

r-^ oo 


Using Theorem 5, (22) and (23) we get the resiilt. 
Theorem 8J If the ^ matrix 


M = 


• 

^l^^n+l-m^ 

• 

p HtI^ 

1 


nf 6^1 

# 


• 

”“P-^o 

m-p-6to 



t: 



is non -singular/ in H^(c^) the coefficients of a 

quadrature optimal multi step forraxila 


(24) 


'n+1* 


^n-s"** ^ ^jr^’'^n+ 1-j' ^n+l-j 


to 


interpolatory for functions approach the 

coefficients b? of the unique multistep formula 


(25) 


yn+l= STn-s* 


j^6. 


to 


2 1 % 

interpolatory for C .. /.../X i 
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Proof: The condition that the matrix M is non-singular 

is sufficient for the forrmola (25) to be unique. ^As both 
the formulae (24) and (25) are interpolatory for » 

in view of Theoren 5 we need to prove that 

lim ^(x^) =0 / i=0/ . . . » 

XT 

where denotes the quadrature error functional for 

( 24) . Proceeding as in Theorem 6 we get 

I Q^(x^) I ^ < Z ^-2(k-i) j ^ 

k=m-pf6^^ 

where is the quadratiire error functional for the 

method ( 25) , Hence 

lim |Q?(x^) 1^ = 0 # 0 < i < 

tr*^ 00 

completing the proof, 

2,7 A Comparison of Different Methods 

In Table 2.8 we present the end point (x^=l-2) errors 
for the set of twenty four differential equations for 
m=l(l)4 step methods of usual and quadrature optimal type 
in H^o^) . Even one step ana implicit methods 

have given better results for equations 1—23, 
equation 24 however, the method has gone worse than the 
usual method. In fact for this equation one step usual 
method has worked better than the other methods of the 
usual type. The 4 step quadrature optimal methods have 
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SIIlMriON: y IN H tc r) ; Rr 2.01 H= .10 
ERRORS AT END POINT FOR EXPLICIT METHODS 
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shown better performance aver the usual ones in all cases. 
The performance of the quadrature optimal methods over the 
corresponding usual methods increased gradually with step 
m=l(l)4. This is as one might expect frcm Figure 2.1 and 
Figure 2. 2 where we have plotted the norm of local trun~ 
cation errors. In Figure 2, l^log^^C | ITJ f ) and log^^C i |T^^ j {) 
have been plotted for the explicit case# and the same is 
done in Figure 2. 2 for implicit case. It is interesting to 
note that the norms of local tiruncation errors for quadrature 
optimal methods ar e all through less than those of the usual 
methods. 

In Table 2.9 we present end point Cx^=l. 2) errors 
for the set of the twenty four differential equations 
dbtained by applying the -methods: expli-cit usual# quadrature 
optimal# quadrature optimal interpolatory for polynonial 
of degree p = 1(1)4 and quadrature optimal interpolatory 
for the functions e*^^ with a = + 1,6, In Table 2.10 
the same information has been given for implicit methods. 

It follows from the Tables 2.9 and 2,10 that the 
quadrature optimal ones have remained the best among the 
usual# quadrature optimal# quadrature optimal interpolatory 
for polynomials (JF# p =1# 2, , . , ) and interpolatory for 
f\anctions exp (l,6x) and e^p (-1.6x) (denoted as M ) -methods 
except in a few Gases(e,g. the explicit method interpolatory 
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for 2 degree polyncinials has performed bettor than optimal 
for equations 13# 14# 15 and 16i the explicit method inrber- 
polatory for exp (l,6x) and exp(-l,6x) has performed better 
on equations 13# 14, 15 and 16 ; the implicit method 
interpolatory for 2 degree polynomials has performed better 
equations 13# 14# 15 and 16 and interpolatory for 
exp(1.6x) and e3pi(“1.6x) has performed better for equation 
10 ). 



C3iAPTER 3 


OPTIMAL MULTISTEP METHODS IN 

3.1 I n-broduct± on 

Quadratxare optimal multistep methods in H^(c^) 
described in Chapter 2 have been seen to be locally inter— 
polatory for fxmctions log(r^-Xj^^^_^x) . In this chapter 
we study the corresponding optimal multi step methods in 
the space H (c^). As will be seen these optimal formulae 
are locally interpol atory for the functions 
the growth of which in the peripheral region of is of a 
higher order than that of the corresponding fxinctions in 
the quadrature optimal case# 

Thus as compared to the quadrature optimal formulae 
of Chapter 2# the formulae obtained in this chapter vdll 
be better suited for differential equations whose solutions 
have more pronounced singxalarities outside c^# 

In Section 3.2 we give a direct proof that the 

2 

derivative of the kernel function of H is the 

complex conjugate of the representer of the derivative 
evaluation functional in H^(cp . In Section 3.3 we 
describe the equaticxis pertaining to optimal mxiltistep 
methods in and thei^r ntmerical implementation. 

Nuraerical results are presented for the set of differential 
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equations described at the end of Chapter 1. In Section 
3»4 optimal methods intesrpolatory for polynomials of 
certain degree are described. In Section 3.5 we consider 
optimal methods interpolatory for a fixed set of functions. 
In Section 3,6 we establish the limiting behaviour of the 
coefficients as r -♦ oo , in Section 3, 7 a ccxnparison of 
different methods is given, 

3*2 The Representer of Derivative Evaluation Functional 

In Chapter 2 we established that the derivative 
evaluation functional at any point inside the disk 
= {z: |z I < r} is a bounded functional in H (c^ , 

As H^{c^) possesses a reproducing kernel function 

(1) K(z,t) = , 

2Tt(r -zt) 

the representer D(trZ^) of the derivative evaluation 
fianctional at a point z^ e ought to satisfy 

= Is K(z ,t) 1^^^^ 

27t(r^- z^t) ^ 

This has been used in Chapter 2 , while calculating 

the expression for the local truncatipn error. For 

> 

the space H^(c^), we directly verify# in the following 
theorem, that D(t#z_) as given above is indeed the 
representer of the derivative evaluation functional# 
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Theorem 1 : The functi on 


( 2 ) 


D(t,2 ) = ^ r SH^CcJ 

° 27t(r^-ZQt)^ ^ 


is the representer of derivative evaluation functional 


at z £ D^, 
o r 


Proof * The only singularity of DCt^z^) is a pole of 
order two at t = r /z , which lies outside C • Hence 

V./ 

2 

the fxincti on belongs to H (c^) , To complete the proof 
we need to show that 


f ‘ iz) = (f (t) , 


rt 


° 2Ti(r^~z t)^ 


By a Cauchy integral formula 




f(t) 


dt 


(t-z^) 


1 

2 %± 


f(t) 


271 

^ ,“ie~772 

o (re -z^) 


re^® ide 


r 

271 


^ f(t) 

,2 - 40,2 

(r -tz^e ) 


r 

271 


/ as 


= (f(t). 


rt 




). 


which proves the result 
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3.3 Optimal Multistep Methods in H^(c ) 


r' 

For determining the coefficients bj of an optimal 
multi step formula 

m m . 

-’"■°to 


where * s are prefixed according to some consistent 
and st-able known usual formula# we have the normal 
equations 

(4) C b = d 


with 


(5) 

and 

( 6 ) 


A ^/\ 

b 


[bg / 


to 


m- 


"•J 


d. = 


[: 


27i;h '^n+l-i ‘ /„2 




m 

- S 


k=l (r^- ^ 


}] . 


i.J = 


which may be solved by Gaussian elimination after 

, /V A 

cancelling the canranon factor r/ 2k in c^ j and d^ « Moreover# 

^ /> /\ 

ronaaoving h from dj^^ results in hb instead of b. 



77 


The following theorem characterizes the optimal 
multistep formula (3) in H=(cp. 

Theorem 2* The optimal multistep formula ( 3 ) in 
in which are prefixed, is characterized by that it 

is locally interpolatory for functions 


y^Cx) 




i = 6^Q(l)m. 


Proof : The proof immediately follows fran Theorem 1«4 
using ( 1) and ( 2) . 

c orollary 1* If 0 some i = 6^^(l)m in (3), 

the formula (3) becomes consistent. 

Pro of: since the corresponding usual formxila is consistent, 

by its precision for constants, S = 1. Further, by 
Theorem 2, if exact for y(x) = x. 

These two facts combined make (3) consistent. 

Analogues of Corollary 1 in the case of quadrature 
formulae have been considered in [30 ] and f 39-40] , 

The norm of the local truncation error 

for a general multi step formula 

m - m 

* ** if* '’J 

a.«i 


is easily seen to be given by 
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C7) IIT^II 


2 r 


n 


[' 


m 

2 a. 


Wn.l-i 


where 


m 


m 


2 a. {-s 3 S a. -R 3 : 

i=l r j = l r 


} 


n-j-l— i 


m 


- h _2_ bj g(x^^^_j) 


j=^ 6 . 


to 


m 


3 ^-' 


-h E S bj ^ 


gCx) = X [• 


m 
S a. 


^ ^ ■ o “"“o 1 • 


In view of Theorem 2# or equivalently the normal 
equations (A), the norm of'T^, the local truncation error 
for optimal formula ( 3 )# simplifies to 


(8) 


A 2 

"V 


r 

27C 




m 

2 

i=l 


i ^2 
r * 


'^n+l^n+l-i 


ro ’ 1 m 3L 

« 2 a. £ —5 ::: ** 2 a^ j 


- h 


J=®to 


In Table 3.1 we present the norm of local trun- 
cation errors and the coefficients of implicit and 
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2 

esiplidt S-step optimal multistep methods in H (c^/ 
where the corresponding usual methods are the Adams - 
Moulton and the Adams-Bashforth. The coefficients have 
the same general behaviour as in the case of quadrature 
optimal methods (i.e., all the coefficients vary fron 
point to point and that, excepti-rq the first coefficient 
of the implicit method/ all the ether coefficients have 
increasing magnitude). 

The local truncation error norm at the beginning 

i.e. at X = “1.5 has reduced (as compared to M case) and 

now it is approximately 1/7000 of the local trucation 

error norm for usual method. The general behaviour of 

NT 1 1 is observed to be as in the case of the quadrature 
n 

optimal methods M of Chapter 2. 

In Tables 3 . Ka) -3. 1(d) we present numerical results 
obtained by using the coefficients given in Table 3.1 on 
the set of twenty four differential equations. We observe 
that the explicit method is just better on differential 
equations 16 and 24 J one decimal place better on equations 
2/8/10,11/12/15/17/18/19/20 and 23; and two decimal places 
better on equations 1/3/4/5,6/7/9/21 and 22, as compared 
with the corresponding usual method. The implicit method 
is just better on equations 15,16 and 17; one decimal 
place better on equations 5, 12,18/19/20/21,22,23 and 24, 
two decimal places better on equations 2,3/4,6/7/8,10 and 11 
and three decimal places better on equations 1 and 9. 
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As compared with the results of Chapter 2^ we observe 
that when the singularities of the solutions are near 
the optimal methods have performed better and in the 
situations where the singularities were comparatively 
far from the quadrature optimal methods have worked 

better, (e.g. for equations 1/2., 3^4/ 5/ 6^ 7,8 and 9 the 
optimal method is better than the quadrature optimal 
method while for equations 13/14/15 and 16 the quadrature 
optimal method has worked better.) 

2 

3,4 Optimal Multi step Methods in H (c^) Interpol atory 
for polynomials 

In Section 1.6 we have desciibed optimal multistep 

methods interpol atory for polynomials of certain degree, 

2 

In this section we describe a H (c^) situation taking 
the points ^ straight line in the domain of 

H^(cp • We asstirae that the points are eq\ii spaced and 

separated by a'dist^ce h. 

Reasoning as in Secticm 2.4 to get an optimal 
multi step foamiula we have to determine the optimal 
coefficients in -the formula interpol atory polyncxidals 
of degree p given by 
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where with prefixed a^*s# s are as in the corresponding 
usual formula. The normal equations are easily seen to he 

(10) = dP 
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In view of Theorem 1,6 we have: 

Theorem 3J The optimal multistep formula (9) is chara- 
cterized by that it is locally interpolatory for functions 
i=l(l)p} u Chj(x); j-p(l)m-6^QL where 
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In Tables 3.2 to 3.6 we present the coefficients 
(calculated via and error norms for H (c^) -optimal 

5-step methods interpolatory for polynomials of degree p# 

1 < p < 5. The general behaviour of the error norms and 
the coefficients has been same as in the case of the 
quadrature optimal multi step methods interpolatory for 
certain degree polynomials, A summary of the numerical 
results obtained by applying these coefficients on the 
twenty four differential equations is as follows: 

As compared with the usual explicit case# the H 
explicit optimal 5-step method with p»l has performed just 
better on equation 24? one decimal place better on equation 
2 , 5, 10 # 11# 12# 18/ 19# 22 and 23? two decimal places better 
on equations 1 # 3 # 4 # 6 # 7 # 8 # 20/21 and 22# and two to three 
places better on equation 9. Similarly# as compared 
with the usual inplicit 5-step method the implicit optimal 
5-step method has performed just better on equations 
15# 23 and 24? one decimal place better on equations 
3 , 5, 12# 18/ 19# 20# 21 and 22; two decimal places better 
on equation 9. In the case of equations 13 / 14# 15# 16 
and 17# however# the performance of both the usual 
and the optimal methods has be^ comparcble. 
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The explicit method vath p=2 has performed just 
better on equations 15,23 and 24 J one decimal place better 
on equations 2,13/15,19,22 and 23; and two decimal places 
better on equations 1,3,4,5,6,7,8,9,10,20 and 21. The 
implicit method has performed just better on equations 
11,12,14 and 24; caae decimal place better on equations 13, 
15,16,17,18,19,20 , 21,22, and 23; two decimal places better 
on equations 1, 2, 3, 4, 5, 6, 7,8 and 10; and three decimal 
places better on equation 9. 

The explicit method with p=3 has worked just better 
on equation 24; one decimal place better on equations 
1, 2,5 , 6,8,9,10,15,16,19 , 20 , 22 and 23? and two decimal 
places better on equations 3 and 4. The implicit method 
has worked just better on equations 5 and 24; one 
decimal place better on equations 5,10,16,19,20,21,22 and 
23| two decimal places better on equations 1,2, 3, 4, 6 and 7j 
and three decimal places better on equations 8 and 9. 

The usual and optimal methods on equations 11,12,13,14,17 
and 18 have performed almost catparably, . 

The ejqslicit method with p=4 has worked just better 
on equations 10,13,14,15 and 24;, one decimal place better 
on equations 2,5,6,8,19,20,22 and 23; and two decimal 
places better on equations 1,3, 4, 7,9 and 21. The implicit 
method has been just better on equations 10,13,14,15 and 24} 
one decimal place better on equations 1,3,4,5,6,7,8,19,20,2,2.' 
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22 and 23;3^'S two decimal places better on equation 9, 
On equations 11 and 12 the performance of usual and 
optimal methods has been almost similar. 


The implicit optimal method with p=5 has been just 
better on equations 10,13/14^15/16,17 and 18? one decimal 
place better on equations 1/2/3/4/5/7/19/20,21/ 22, 23 and 
24/ two decimal places better on equations 6,8 and 9, 

The optimal and usual methods have worked almost similarly 
on equations 11 and 12, 

2 

optimal Mvatistep Methods in H (c^.) I nterpol atory 
for a Set of P reassigned Functions 

The coefficients "bj of an optimal multi step formula 

m p , x 

^n+l* , ^-®i^n+l-i* ^^^+l-j'^n+l-j^ 

J = \o 


3.5 


( 11 ) 


in H^(c ) with prefixed a^'s interpolatory for linearly 
independent functions £9^/*../^^} are to satisfy the 
normal equations 
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with C/'d as in equations (4), \ 
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(14) 
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and 0 is a p X p nvill matrix. 

This method is characterized as follows: 

2 

Theorem 4: An optimal multistep method in H. (c^) inter- 
polatory for linearly independent f^lnctions ^ (Pj,/ . . • ^ 
is characterized by that it is locally interpolatoxy for 
functions 

U {hj^Cx)? k =6^Q(l)m-p } 

with 


hj^(x) = 


m 

+ 2 


X 


/ 2 - ^2 . ir Z -sr't 

"^n+l-k^^ x=sm-p^l ^ to (r 


where G 




=l(l)p 


“1 

with G = W E , 


q=l(l)m-p+6^^ 


where 

W ss C^i i=l(l)p 

j=m+l-p(l)m 


and 

E = C''’h=si»i-j>^i=i(i)p 

Proof: The proof follows from Theoraii 1*8. 

In Table 3.7 we have displayed the optimal case 
interpolatory for functions '^(x) = exp(1.6x) and 
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92,(x) = (“1.6x). The behaviour of error rtorms ahd that 
of the coefficdents is found to be as in other situations* 
The performance of the methods as compared to usual methods 
is as follows: 

The explicit optimal method in H^(c^) has been 
just better on equations 14,17,18 and 24; one decimal 
place better on equations 2 , 10 , 11 , 12 , 13 , 15 , 16 , 19 , 20 , 21,22 
and 235 and two decimal places better on equations 1,3, 4, 5, 
6,7,8 and 9, The implicit optimal method was observed 
just better on equations 14,18,23 and 24* one decimal place 
•better on equations 12,13,15,16,17,29,20,21 and 22; two 
decimal places better on equations 1,3,4,5,6,7,8,10, and 
11; and three decimal places better on equations 2 and 9* 


3* 6 Behaviour of the Coefficients as r-* 

In Section 2.6 we have studied the properties of 
the coefficients of quadrature optimal multistep methods 
in H^(cp. Here we present the corresponding results for 
optimal multistep methods in H^(c^) • We first prove: 


Theorem 5: Let 


(15) 


m 

,« 2 a,y, 

i=l 


^n+1" j “i^nfl-i" 


^ F \ 


he a mill ti step formiola interpolatory for functions 

^95 and let the matrix [‘^i^^n^l-j^^i, j«6_ 

to ^ 

be non -singular. Let 
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( 16 ) 


m ^ F / 

i'=l J“%o 


.) 


be an optimal mxxLtistep formula xn H (c^ interpola ry 

.CD 9 ^ }, 0 <P <m-6to 

for . If the local 

non-lnteSolatory case for tf^ese fmctxons) . I* th 


truncation error 

(17) 


functional tJ, for C16) satisfies 


lim ° 3 = P*®toU> 

XT 


m 


then 

( 18 ) 


and 


/V 

lim b 


S' = b? , i = Sto^""- 


r-* °® 

Proof 1 Let 


Jt 




FnT 


B' 


'» [Pa 


"to 


« A p nT 

A i> fO* F .b'*^ t 

■Q^sz Lbe r **** ' mr-' 

t 6to^ 


T, 


nr to 

^en, substituting tbe actions 

^, 3 (15) ana (16) subtracting 

the formxxLae ^xo; 


we get 


hG •<• V “ ° 

^ =0 by (17) and as 

Since li^ ^nr ^ 

result is immediate. 


G -is non-singul ar^ the 



Theorem 6: Let a. ^ i=:l(l)m/ be constants such that 

E a.= 1. Then the H (c ) -coefficients b. of an 
x-1 r jr 

ox:)timal multi step formula 


(19) 


m 


m 

2 


A 


~ -^“^to 


approach the coefficients b^ of multistep formula 


( 20 ) 


m 


m 




j=6. 


to 


of maximal polynomial precision. 


Proof : Since the formula (20) is interpolatory for 

polynomials of maximal degree and it is assumed that a^^'s 
are known constants with 2a^=l/ we know that the coeffi- 
cients b. of method (20) can ■uniquely be determined by 

■J 2 

using that it is interpolatory for functions {x,x ,...,x 

In view of Theorem 5 we need to prove that the truncation 

error fiinctional T^^^ for (19) satisfies 

lim T (x^) = 0 , 0 < i < m+6^^ 
r*»«» 

The method (19) being optimal, we have 

II V" i "L" 

where T is the local truncation error functional for 
n 

(20). Since SS]^Co' as defl-nedin (2.3), is a complete 

2 

orthonormal sequence in H (c^,), we have 
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i=0 ^ k=0 ° ^ k=ro*6^,.n2^ " 


Thus/ for each fixed i^ 0^i^n\ + have 


yx 4 O «0 

1 T (x^)l ^ < S 

k=m+6^^+l 


^-2(k-i) 1 2 


The e3(pression on the right approaches 0 as «> . Hence 


lim = 0 / 0<i<m-f 6^^ ^ 


JC-* oo 


tl 


which proves the result, 

m 

Theorem 7, Let S a. = 1 where a. 's are given constants, 

i=l ^ ^ 

2 P 

The H (c ) -coefficients b. of an optimal mtiltistep 

X» J 

formula 

( 21 ) 


m 

2 

iftl 


m 


<LP 




j^6. 


to 


Interpol atory for polynomials of degree p <m + 6^^, 

approach the coefficients of the corresponding usual 

formula as r • 

m 

Proof: Since S = 1 the formula (21) is interpolatory 

i^tl ^ 

for constants. As it is interpolatory for polyncxnials of 
degree p# the truncation error functional of (21) 

satisfies 

"t?^ (x^) s 0 ♦ 0 < i ^ P , 
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In view of Theorem 5 it suffices to prove that 


lira ~ ^ ^ ®tl 

IT-^oo 


AP 

I imXr 


Since I IT^ J1 < 11 I! where T„ denotes the truncation 
nr n n 

error fionctional of a usual method/ as in Theorem 6, 


we have 


2 r-^lT^ (x^) < S r'^lT^Cx’^) 

i=p+l kiST\+6^^i-l 


Hence 


. <i,P / i s , 2 „ 

lT„„(x ) I < S r 

ks3n+6^jl^+l 


-2Ck-i) j2^ p+1^^+6 


tl 


Since the right hand side approaches 0 as r-*'® we have 
lim =0 / P +1 < i < m +6 ^ » . 


nr 

00 

completing the proof. 

m 


Theorem 8 i Let S a^,. =1 where a^'s are given constants, 
ial ^ 

Let the (^ 1 + 6 ^ 3 ^) x (m+ 6 ^ 3 ^) matrix 


M 




1 


"tl 


m-p~ 6 , 


^+ & 


to 


tl 


'^i^^n+l-m^ 


^p^^n+l-m^ 

1 

^+l-m 


^•♦•1-m 
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be non-singular. Then, the H (c^) -coefficient of 

an optimal multistep formula 


( 2 . 2 ) 


^ A p 

■^“°to 


interpolatory for functions p < m + 6^^ 

approach the coefficients bj of the tmique multistep 
formula 


(23) 


m 


m 




to 


2 t 

interpolatory for functions C 9^, . , . ,x * 

Proof : The condition that the matrix M is non-singular 

is sufficient for the formula (23) to be determined 
uniquely. Since both the formulae are interpolatory for 
functions f9-,..,,9 3 , in view of Theorem 5 we need to 

ir 

prove only that the local truncation error functional 

T^^ for (22) satisfies 
nr 

lim T^ (x^) =0, 1 < i < m-p-*-6^j^ . 

r-* «> 

A similar procedure as in Theorem 6 gives 


m 


IT^ (x^)l^ < 2 r’"^^^“^^lT^(x^)l^, lo<m-pi*6^^ 

” ]c=!m-p+ 


The right side of the above inequality approaches zero 


as r •* . Hence 

/-i 


lim t;:^(x ) *= 0 , 1 < i < m-p+6. , 

HIT 


i.ihi ch nroves the result* 
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O 

3,7 A Comparison of Optimal Methods in H (c^ 

In Table 3#8 we present the end point 2) errors 

for the set of twenty four differential equations for 1 to 
4 step e 3 <plicit and implicit methods alongwith the 1 to 4 

f 

stop usual methods. We observe that the performance of ^ 
the optimal methods has been better all through except j 
in the case of equation 24/ which seems to be a chance. 

/s 

In Figure 3.1 we have plotted log^^C I !T^l 1 ) and 
log,,.( II T 11) at different points for explicit methods of 
step 1 to 5, We observe that the error norm for the 
optimal method/ as expected# remains throughout less than 
that of the usual method. In Figure 3.2 these norms have 
boon plotted for the implicit methods. 

In Table 3.9 we present the end point errors (x^-1. 2) 
for the twenty four differential equations for 5 step 
explicit H^(Cj.) optimal method/interpolatory for polyno- 
mials of degree 1-4 (as denoted by the 

superscript digit denoting the degree of the polynondals 
for which it is interpolatory)/and interpolatory for 
e 3 ^( 1.6 x>/ expi-U6x} (denoted as 

We observe that has performed best in many 
cases (e.g. explicit method for equations 9.13/14/15/16/ 
19/18 and 23 ) andft^cm other equations (e.g. for equations 
2 / 3 / 4 / 5/6 and 7 ). ^er optimal methods have also performed 

best on selected equations# 



Table 3.8 


SITlUnON: y IN H (c r) J Ra 2.01 Hr .10 

errors at end point for explicit HETHDdS 


Eq.No, 


1 ? 

il 

16 

17 

18 
19 

21 


24 


USUAL 

zBlBl 

0.12E+00 

-0.44E-02 

•0.50E-03 

-O.eiE-02 

-0.20E-02 

'0.17E+01 

-0.13E+g4 

0.13E+54 

-O.99E-02 

-0.85E-02 

-0.80E-02 

-0.69E-02 

-0.67E-05 

-0.66E-54 

0.63E-01 

-0.60E-02 

■S:55l:8? 


^OPTIMAL 

-0,25E-0r 

-o,l9g;-o| 

zkmli 

0,66E-01 

-6,30E-02 

-0.36E-03 

-0.55e-02 

-0.14E-02 

-O.lSE-fOl 

■8:88l:8| 

:8:Jp:8| 

-0.43E-02 

-0.36E-02 

0.43E-02 

-O.lOE-01 

-O.llE-01 

“0.46E-02 

-8:Ui!8§ 


_USUAL_ 

0.2ic-02' 

8:l!|:8i 
.8;il|:85 
8:l?|:8i 
8:l5|:8i 

0.i7E+0Q 

0.76E+03 

■8:!?|:8I 

':K|:8! : 


0.16E-03 

0.82E-06 

0.79E-05 

0.53E-02 

8:5P:8i 


OPTIMAL 

0.24E-g4' 

•0.69E-04 

0.26E-03 

0.57E-64 

0.20E-03 

0.14E-04 

0.52E-05 

-0.55E-04 

0.89E-06 

0.58E-03 

O.2eE+00 

-o.eeE+oi 

-0.34E-04 

0.26E-04 

'0.40E-04 

•0.33E-04 

-0.15E-06 

-0.16E-05 


^osual 

0.3QE-02 
‘47E-g 
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In Table 3.10 we present the end point errors of 
5~stGp implicit methods. Here we have a similar obser- 
vation as in the case of explicit methods. Thus optimal 
methods intexpolatory for polynoraials/scxne other functions 
have performed better than optimal methods in selected 
situations. 



CHAPTER 4 


QUADRATURE OPTIMAL MULTI STEP METHODS IN L^(cp 


4,1 I nt roducti on 

2 

The noOT in the space H (c^) is based on a line 
integral along the boundary which takes care of the 
points in the interior of c^ only indirectly* Hence 
the points in the neighbourhood of the boundary receive a 
dominant weightage in the mrdtistep methods corresponding 
to the space H^(cp. As the maxim\jra use of these methods 
is envisaged on a subinterval of (-r,r) it effectively 
means that the H^(c^) methods would# relatively# be 
more effective in the peripheral regions. On the contrary# 
in the case of the Hilbert space L^(c^)# however# the 
norm corresponds to an area integral on the disk and 
thus it intrinsically gives relatively more weightage to 
the points in the interior of c^. The corresponding 
multistep methods are conseguently ej^ected to result in 
a better distributed behaviour in the central regions. 

It may also be noted that (as sets) L (c^.) H (c^. 

This ch^ter is a study of quadrature optimal 

T 2/C \ Tn Section 4.2 
multistep methods in the space L \c ^) . 

^ L^(^ ) and establish the boundedness 

we define the space n ^ ^ 

of derivative evaluation functionals in L (c^). m 



Ill 


Section 4.3 we present the quadrature optimal multistep 
2 ^ 

methods in L (o^) . in Section 4.4 quadrature optimal 
multistep methods in L^(c^), interpolatoary for polynomial 
of a partictilar degree# are presented. In Section 4,5 
wo describe quadrature optimal m -step methods interpol atory 
for arbitrary functions. In Section 4.6 we discuss the 
limiting behaviour of the coefficients and in Section 4,7 
we present a conparison of various quadrature optimal 
multistep methods in L^(c^). To distinguish from the 
H^(c^) situation# in the notations we have used an under 
bar (c.g,# as in b^), 

2 A 

4.2 The Space L (c ) and Boundedness of Derivative 

XT 

Evaluations 

2 

Let L (D^) denote the space of complex valued 
measurable functions h{x#y) defined on {(x#y) : 

<, for which Ihl is square integrable in D^, 

The space L^CD^ turns cut to be a conplex Hilbert space# 
with the inner product 

(h,X) « // h(x#y) k(x#y) dxdy. 

D 

r 

The subspace of L^(D^) consisting of functions 
of the type hCx,y) = f(z) (z = x-i-iy)/ vAxere f(z) is an 
analytic function of z# regular in D^# rs closed and has 
an infinite, dimension. This is the Hilbert space L ic^) 
with the induced inner product given by 
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=// f(s) g(z) dxdy • 

D 

r 

The space L (c^) has a reproducing kernel 
function as 

2 

(2) K(z,t) = £ 

TtCr -2t) ^ 

and is separable with a complete orthonormal sequence 
of functions given by 

k 

(3) (s) = - , k = 0,1,... 


We prove that derivative evaluation of any order 
at a point inside is a bounded functional: 

Theorem 1: For z e and k any positive integer 

is a bounded linear functional in L ('c^) with 
2 ^ 

k+ — 

. . _k. , . . , . , r”' / 2 2 

1 I D (z)|J < klNT^ ( 


Proof . By a cauchy integral formula we have 




, Izl < p < r. 


, .k-^-l 

Cp (w-z) 


Let 


|z 1 » r-6 and ^ = Pe 


i© 


Then 




Integrating this inequality with respect to p from r-69 
to r, (0 < ? < 1) ' we get 
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6<P <|i / ,“jfi£^ papde 

p=r-69 9=0 (9“r+6)^'*'^ 




Stc 


PdPde 


1/2 


p-r-69 0=0 (p~r+6) 


2k-i‘2 


If 


t O 

L (cp 


Now let us evaluate the integral under the square root. 

r 2% 

/ — -P.dPde 


p=r-<59 0=0 (p-r+6) 


= 2ji / 




p=r-6q> (p-r+6)^’^2 


2t [ 


r - 69 


2ki*l {(1-9)6 }2ktl 2k+l ^2^+1 -» 

r 1 ' 

^ ^ ( 2k^l) (2ky ^ (p-r+6) 


(2k-»-l)6 


^2-^ {(r_69)(l_<p) C2k+l)^^j 


2n 


[ l] 


( 2k+l) { 2Ic) 6 


2k 


2iir 




Choosing 9 » vre get 




i W'f|f 
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Hence 


I Id1!| I < 

2 ■“ 2 % V-I-t: H / 


271 ^r~tz 1' 
which proves the result, 

4,3 Quadrature Optimal Multistep Methods in L^(c ) 


The coefficients bj in a quadrature optimal 


2 

formula/ in L 

( cp . 



satisfy the normal equations 

(5) 

to ro 

C b = d 

with 



r2 

%J 


and 


4 

- ^-s 

f r^- X X ) (r^ 

^n+1 ^+l-i^^^ 


i#j * 

The following theorem characterizes such a multi step 
founxila. 

Theorem 2, The quadrature optimal multi step formula (4) 
is characterized by that it is locally interpolatory for 


functions 
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Yi (x) = 


> if ^., 1-1 ^ 0 


X , if X 


n+l-i 


0 . 


^ ~ "^to' 2, 

Proof, Proceeding as in the proof of Theoron 2. 2 and 
malting use of (2) wa find that the foimula (4) is inter- 
polatory for the functions 


yi(x) = 


fl rj 1 

i " ^+1-1 


' i* ^,1-1 ° 


I - 




to' 


from which the result follows. 

It is worth noting that if r,h and in the 

normal equations are respectively replaced by ar, oh and 
a normal equations do not change and hence 

the coefficients rfsnain the same. This implies that 

working with one value of r and calculating the coefficients 
b . for different values of h is sufficient to give rise 

''J 

to the coefficients for different values of r. This 
observation is easily seen to be applicable to the methods 
of the previous chapters as well. 

The L^(0j^)-norm of local truncation error functional 
^ at point x^^i for any multistep formula of quadrature 
typo with coe££ici6nt$ bj turns out to be given y 
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(6) nT„n2 = f .2,-2 


n-*s' 


2 ^ 

-2Ref(r j; b g(^ ) 

to 


m 


to 


-2h b. ,] 

^“^to *'^+l-i^+l-i^ 


where 


g(x) = 2x - (r^-x^_gX)“^} . 


-3, 


In the numerical results to follow (5) and (6) have 
been implanented for L^(cp with r = 2,01 and corres- 
ponding to the standard situation of h=,l with 3^= ih/ 
i=0# +1#*,. of earlier chapters. 

Table 4,1 displays the local truncation error norms 

and the coefficients of L^Cc'^) -quadrature-Optiraali (mtS) 

multistep formulae in the implicit and explicit cases. The 

general behavioxir of coefficients has been the same as in 

the case of H^(cp formxalae. The range of variation of 

coefficients in the present case has however , slightly 

increased. At the point x = -1. 5 the error norm in the 

r L ^(Cj^) —quadrature optimal case is about 1/7500 of the usual 

I case. This indicates that the norm of local truncation 
I 2 

' error functional has reduced as compared to the H (c^)- 

1 quadrature optimal case. 
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Tables 4. 1(a) ~4, 1(d) are applications of Table 4,1 

on the set of twenty foxir differential equations. The 
2 

^ quadrature optimal 5— step method has 

performed just better on equations 14,16 and 24; one 
decimal place better on equations 2,4,6,8,11,12,13 ,15, 
17,18,19,20, and 23 and two decimal places better on 
equations 1,3,5,7,9,10,21 and 22 as compared to the usual 
method. The implicit quadrature optimal method has been 
one decimal place better on equations 3,11,12,13,14,15,16, 
17,18, 19, 20,21, 22, 23 and 24; two decimal places better 
on equations 1, 4, 5, 6, 7,8,9 and 10; and three decimal 
places better on equation 2, as compared vd.th the corres^ 
ponding usual implicit method. 


2 

4,4 Quadrature Optimeil Multi step Methods in L‘^(c^) 
Interpolatory "for Polynomials’ 

As in Section 2,4, for equidistant points 
•~-P 2 ^ 

the coefficients Y- of a L (c ) -quadrature optimal mxiltx- 

“TJ ^ 

step f ormul a 


(7) 


7’n+l“^n-s " jlo ^ 


P-1 u -i 

-h E r, ) 


nt 


\l 


"tl 


m-6 


to 


= h S f(x 

j.p ^^-tl 




interpolatory for polynomials of degree p < m -<■ are 

sxitiilarly to satisfy the noannal equations 
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( 8 ) 


c r 




where / for i# j = p(l)ni-6+.^' 


and 


I - ~ '■ Z - n '* ^ » 


to 

T 


—p' * * * 


to 


I ” S (-1)^*^(^)(^) 
^ 1^0 ^ ^ 


d? = ~ [ 2 (-l)^(^) — ^ 


/ 2 
(r -X 


‘n+ 6^ -k^n+6^ ^ -1 ^ 


^n+1 ~ ^n-s 


"i nh 


1=0 




tl' 


■'tl 


"'^2%" cl 2 (-l)^^^(i)(2) 

^ Cf-l y-V-l_l-N XJV 


a:i=0 ^ 'k=0 1=0 




The following theorem characterizes fonn\ila (7) • 

Theorem 3 : The multi step formula (7) is a guadratiure 
optimal formula if, and only if, it is locally interpola- 
tory for the functions 


J k=0 


}* . 


hj(x) = 2 (-1) CjJ.) i r ,2 - 3 jX 

j = p(l)m-6^Q, 


2 . 


= 0 . 


where { , }* is to be replaced by x/r if 

The proof which uses (2) and is similar to that of 
Theorem 2.3 is omitted. 
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The multi step coefficients bj in the following 
numemcal results are obtained frcm the coefficients 
j as determined through (8) . 

In Tables 4,2 to 4.6 the coefficients and error 
norms of L ( c^) -quadrature optimal 5-step methods Inter- 
pol atory for polynomials of degree 1 P < 5 are displayed, 
A summary of numerical results obtained by applying these 
methods on the set of twenty four differential equations 
is as follows! 

As compared with the usual explicit casesthe L^(c^)~ 
quadrature optimal 5-step method with p=l has worked 
just better on equations 13# 14 and 24? one decimal place 
better on equations 2,6#8#10# 11# 12# 15# 16# 17# 18# 19 and 23# 
and two decimal places better on equations 1#3#4#5#7#9, 

20# 21 and 22, In the implicit case the method has worked 
just better on equations 14 and 24# one decimal place 
better on equations 3# 11# 12# l3# 15# 16# 17# 18# 19# 20# 21, 22 
and 23 > and two decimal places better on equation 2, 

The e 3 q>licit method with p=2 performed just better 
on equations 13# 14# 16 and 24# one decimal place better on 
equations 2# 6#8# 10# 11# 12# 15# 17# 18# 19 and 23? and two 
decimal places better on equations 1# 3#4,5#7#9# 20# 21 and 
22, The implicit method has been just better on equations 
12/18 and 24 J one decimal place better on equations 
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SITUATION! Y IN L, (C r) ; R= 2.01 H= ,10 

HULTI STEP METHOD OF 5 STEPS INTEftPQLATRf UpTO 3 DEGREE POLYNOMIALS 
"eOPTIMAL BU) H X BC2) H X B(3) 
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11, 12, 13# 14, 15^ 16, 19# 20# 21# 22 and 23jand two decimal 
places better on equations 1, 2, 3, 4,5,6, 7,8,9 and 10. 

The explicit method with p=3 has been just better 
on equations 11,12,13,14,16 and 245 one decimal place 
better on equations 1,2,3,5,6,8,9,10,15,19,20 ,22 and 235 
and two decimal places better on equations 4,7, and 21, 

The implicit method with p=3 has been just better on 
equations 11,12 and 24J one decimal place better on 
equations 3,10,13 ,14,15,16,19, 20,21, 22 and 23 and two 
decimal places better on equations 1, 2, 4, 5, 6, 7,8 and 9. 

The two methods have been comparable on equations 17 and 18 

The explicit method with p=4 has been just better 
on equations 14,17,18 and 24* one decimal place better on 
equations 1, 2,3,6,7,8,9,10,13,15,16,19,20,22 and 23 and 
2 decimal places better on equation 9. The implicit 
method with p=4 has been just better on equations 15,17, 

18 and 24i one decimal place better on equations 1, 2,3,4, 
5,6,7,8,10,13,14,16,19,20,21,22 and 23? and two decimal 
places better on equation 9. The methods have been just 
worse on equations 11 and 12, 

The implicit method with p=5 has been just better 
on equations 10,13,15,16,17 and 18; one decimal place 
better on equations 1,2,3,5,6,7,8,14,19 , 20 , 21,22,23 and 

better on equations 4 and 9. 


24 and two decimal places 
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4.5 Quadrature Optimal Multistep Methods in 

Intexrp olatory for a Set of Preassipned Fvinctions 

The coefficients Bj of a quadrature optimal multi- 
step formula in L^('cp 


(9) 


m 




interpolatory for linearly independent f\anctions 
{9-,. ..,9 }* are to satisfy normal equations 


CO 

C 





1 

1“'°' 


0 




e 

mssA 


where C, d are as in equations (5), b =Lk * 


] t 


to 




and 0 is p X p null matrix. Here we have: 

Theorem 4; A quadrature optimal multi step method in 
L^(c ) interpolatory for functions {9- #...,9 } is 
characterized by that it is locally interpolatory for 
constants and functions 

{93^... ./ 9 p} U {hj^Cx), k = p+S^^CDm} , 

where 
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h„(x) “ 


P-ll-6to P-l-"6to 

+ s f 2 w . 

J *6 






to 




P ^T, -l-j 1 ' 


-1 


where 


W = 


r w 1 

•- ij-*i=l(l)p 


to 


with 


W 


-1 


I- ’’h^n-i-l-jhj^KDp 




and -it \jnder tho bracket means that the esqjression witliin 
is to be replaced by ^ if 1=^ l^j . 

Proof , The proof is immediate from Theorem 1,3 using (2) 
after simplifications. 

In Table 4,7 we display the local truncation error 
noOTS for L 2 “quadrature optimal multistep explicit 
and implicit methods interpolatory for functions 
9j^(x) = exp(1.6x) and ^ 2 ^^^ “ exp(-l,6x) and the coeffi- 
cients of the methods. The general behaviour of the 
coefficients has been as earliec. The nxomerical resvilts 
obtained by applying the coefficients of Table 2,7 on the 
set of twenty four differential equations are summarized 


as follows: 
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As compared to e3<plicit usual method the explicit 
quadrature optimal method has performed just better on 
equations 13/14/15 and 24J one decimal place better on 
equations 2,3/6,8,9,10,11^12,16,17,18,19 and 23; and tv/o 
decimal places better on equations 1,4,5,7,20,21 and 22, 
The implicit quadratixre optimal method has performed just 
better on equations 11,12,17,18 and 24; one decimal place 
better on equations 3,19, 20, 21,22 and 23; and two decimal 
places better on equations 1,2,4,5,6,7,8,9,10,13,14,15 
and 16/ in comparison with the corresponding usual implicit 
method. 

4,6 Limiting Behaviour of the Coefficients as r-*oo 

In Section 2,6 we studied a limiting behaviour, 

as r oo , of the coefficients in quadrature optimal 

2 

multi step formula for the space H (c^) • Analogues of 
the results obtained therein also remain true for the 
space L'^(c^) with little change in the proofs. 

Accordingly, in this section we merely include the state- 
ments of the L^(c^) res\alts and emit their proofs, 

the only essential ^fference being that instead of the 

2 

orthonormal basis defined in (2,3) for H (c^) , 

one is now to use the basis { ^ ^*"r^ defined 


in (3) of this chapter. 
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Th eorem 5 : Let 


m 

3=5 




to 




-oo a quadrature type usual mill ti step formula Interpol atory 
for functions ^ ®5to^ * • ^ matrix 

[ \ i-6 be non -sing alar. Let 


( 10 ) 


^n+1 


= ^n-s-^ ^ 


J"5to 


“jr 


be the quadrature optimal multi step f omul a in L'^(c^) 

interpol atory for functions ,,,.,9 3/ 

to . p-ito^^ 

5. P S meaning a non -interpol atory case). 

If the quadrature error 2^^ for (10) satisfies 


lim ^ P’**^to^^^”' ' 


then 

lim bT = b? / j = 6^Q(l)ni. 
r-* 00 “* 

Theorem 6! In L^(c^), the coefficients Ej of a 
quadrature optimal multi step formula (4) approach the 
coefficients of the corresponding usual formula as r ~ , 

Theorem 7 j In L^(c^)/ the coelficionts bj^of a quadrature 


optimal multi step formula 



m 


~P 


^ns*! ~ ^n— s ^ 


•^~ to 

interpolatory for polynomials of degree p < ' 

approach the coefficients of the corresponding usual 
formula as r ^ oo , 

Theorem 8: If the (m+6^^) x (m + matrix 


M = 




tl 


CP' (x 


) 


X. 


’nHrd 


(x ) 


9' (x ) 


X. 


'tl 


h+l-m 




m~p-6 
n’«*6^ 


to 


m-p— 6 


to 


"TO 


is non~singular, then, as r , the coefficients 


b? of the qaaadrature optimal multi step formxila 

«3 ^ 

m ^ p 

yntl = ^n-s* .f . ' 

to 

in L^(cp, interpolatory for fxmctions {9^(x) . . .9p(x) 3 , 
approach the coef£Lcients of the corresponding unique 
multi step formula 

in 

■J” to 


f 

interpolatory for ^ 9^(x) , , . ./^p(x) ,x,x 


i m m m i 


m 

• #x 


-p+6^1 
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-QQ^^Pari son of Different Methods 

In Table 4.8 we present end point (x^=l. 2 ) results 
obtained for quadrature optimal methods in L^(c^) of one 
v-o four steps and compare these vrith the corresponding usual 
, methods. VJe observe that the quadrature optimal methods 
have invariably remained better than the usual methods 
except in the case of equation 24. 

In Figure 4,1 we have plotted log^^ ( I ! T^(x) II) and 
iP^io^ I I ) for the explicit methods. In Figure 4.2 

the same is done for the implicit situation. We observe 

that the general behaviour of the errors is same as in the j 

2 ' 
case of H (c^ situation of Chapter 2. 

In Table 4,9 we display end point (x^=l,2) errors 
for the set of twenty four differential equations corres- 
ponding to usual# quadrature optimal M# quadrature optimal 
methods interpol atory for polynonials of degrees 1 to 4 
(denoted as ,,, and quadrature optimal method 

<vp 

interpol atory for exp (t l, 6 x) (denoted as M )for the 
explicit case. We observe that tne errors are least for 
M on equations 3 # 4 # 6 / 8 # 19 and 22 ; Sor on equations 6,11 
and 12; for on 1,2,3,4,6,10,15,16,17 and I 85 for on 
equation 237 for m'* on equations i and 24; and for M on 
equations 13,14,20 and 21. From Table 4.10, where we present 
the corresponding implicit situation, we observe that the 



Table 4.8 



, 2 * 

SITUATION! Y IN t, (C r) ; R» 2, OX H= .10 
ERR0RS_AI^END_P0INT_F0H_EXPL1CIT METHODS 

- - — _ 

Eq, NO. H 

USUAL OPpMAL ^USUAL OPTIMAL USUAL OPTIMAL U^UAL OPTIMAL 

“*“7 -0.41E-01 -0,32E-oI 0.25E-02 0 .26E-03*-0r30E-02'"-ol87E-03**0ri / e7o2 

:8;io|:8! 8;i§p! 8:S2|:8S :8;4’4|:?1 
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Figure 4.2: Local Toruncation Error in Implicit L (c^) 
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errors are least for M on equations 7*11/12,18/19,20/ 

««1 *^2 

21 and 24; for on equations 6 and 18, for M on 

equations 10 and 17; for on equations 1,4,5 and 8, 

for on equation 23, for on equations 4 and 22, 

and for on equations 2/9/13,14,15 and 16, 



CHAPTER 5 


OPTIMAL multistep METHODS 


IN L^(c ) 
r 


5,1 Introduction 


In this chapter we con^-i^ 

e consider optimal multi step 

formulae for the Hilbert sn;.r^« t 2.-'''- ^ 

pace L (c^) of functions. 

On the whole, the multi 

^xuLstep formulae of the present 

chapter have the same 

advantages over the corresponding 

quadrature optimal formulae in t 2/1 n ^ 

ae in L (c^) as the formulae 

of Chapter 3 have over those o-f ^ ^ 

wiuse or Chapter 2. Concerning 

the relative growth on (--■r r'^ ^ 

^ V r,r; of functions to which the 

formulae are to be applied, a look at the respective 
interpolatory funoUons reveals that the formulas of the 
present chapter admit functions of a higher growth than 
those of any of the Chapters 2-4. ihls is, of course, 

as might have been ejpeoted theoretically from the inclu- 
sion H^Cc..) L^(c ), 


In this chapter in Section 5.2 we directly verify 

that the representer of the derivative evaluation 

2 A 

functional in L (c^) is the complex conjugate of the 
derivative of the kernel funcUon of L^(cp. Section 
5,3 contains normal equations for constructing optimal 
multistep formulae in a characterization of the 

formulae and numerical tables and results obtained for the 
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• 2 ^ 

same. In Section 5,4,L (c^) -optimal methods interpolatorv 
for polynomials of a certain degree are presented. In 
Section 5.5 we consider L (c^) optimal methods interpola- 
tory for a set of prefixed functions. In Section 5.6 
limiting behaviour of the coefficients of L ^ -optimal 
fonaulae/ as r -» oo / is described. In the last section 
different methods of this chapter have been compared for 
their numerical performance. 


5,2 The Representer of Derivative Evaluation Functional 


In Chapter 4 we proved that the derivative evaluation 
functional at a point inside c^ is a bounded functional 
L^(c^). In the following theorem we obtain the representer 
of the same. We recall here the expression for the kernel 
fmction of L (c^) J 


( 1 ) 

Th oorom 

point 2 


K(2/t) » 5 — zr2 * 

7i(r -2t) 

1: In the space L^Cc^ / the representer 

of the derivative evaluation functional at 

inside is given by 
o ^ 


a 


( 2 ) 


Proof. 


D(t,2o) 


2r\ 

n(r^-tz^)^ 


2 A 

LetfeL (c^. 


Then/ for a inside c^ , 


we have 
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(f(z), ) 


ee fiz) 2 r^z . _ 

■* J ^— 2 — -r^ dxdy 

D„ % ir-zz_) 


o 


r 2 % 10, ^ 2 -18 

^ pdpde 

p=0 0=0 7t(r^-pe 


" i ? U£pJ 2 r^ . (£^li 

It(£5Pe«-z)3 P 
pZ o 


p=0 0=0 


= r { f 


fisiz 


P =0 c . , V 

r %x{z ) 


rhie^/r^2l ^ i ap 

3 p^ 


r 

f 

p =0 


2 p^ 


/ rf(z)zT' 2 ^ / 

p _ r^ 

2 = o 2 ^ 


by a Cauchy integral formula (for second derivative). 
If 2 ^= 0 # the last integral equals 

XT 3 

/ [ f" ( 0 )* 0 + 2 f ( 0 ) ] ^ dp = f’ ( 0 ) . 

p =0 r 

On the other hand# if 0 # the integral equals 

Zq 

/ [£( 2 )z]‘' dz 

z =0 Zq 

z 

= ^ [iz[f(z)z]'JL!o [f(z)z3‘dz] 

Zq 

*» jE * ( * 
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Thus the proof is complete. 

Note that (2) la equiva lent to saying that 
N{t,z ) = 0K{s,t)/3z)l 

2=Z 

O 


5,3 Optimal Multi step Methods in L^C'cp 

The coefficients of an optimal moil ti step formula 


(3) 


m m 

^tO 


with a .^ ' s chosen according to same hnown usual method, 
satisfy the equations 

. /'v A /S 

(4) C b = d 

v/ith b = [ b * **‘*'-n I ' 

""^to 

(5) , 

%j 71 ,2 - ,4* 

and 

]c=l 

The following theorian characterizes the multistep foiroula (3). 

Theory 2: The optimal multistep formula (3) is locally 

interpolatory for functions {x(r ' i=6^Q(l)m}, 

Proof. The proof is an easy consequence of (4) or of 


Theorem 1,4 and ( 2) 
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It may be noted, that if x = 0 fon some 

n-i-l — 1 

^ ^ then the formula (3) is exact 

for polynomials of degree one. 


We also observe from Theorem 3 that the growth of 

the interpolatory functions for (3) is according to the 

third power of The same for quadrature 

2 o 

optimal H (c^) t optimal H'^(c.^) and quadrature optimal 
L (c^) is according to log (r “Xn+l-i^^ ' 

2 -I 

and (r / respectively. This substantiates 

the remark made in the first section of this chapter. 

The L^Cc^)-norm of the local truncation error T^ 
at ^^2. ^ general multi step formula 


m 


is:n+l-i J-'Sto 

is given by 

(7) I lEnl l" = r 


^ 2 ^ ry ^ 


m 


- S 5.{(r 

x^l 


-2 


^ o ^ — 2 , / 


m 




j=6 


to 


m 

-h s 

ks=6to 

m r^+2x^+l-i^n+l--k 

-2h S b. — ^ -‘r ^ 

j-\o "^+l-j^n+l-k^ 


J] 


# 
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where 


( 8 ) 


9(K) = 


Using the noimal equations (4), the noim of the 
local truncation error ^ for the optimal formula (3) 
con be simplified as 


(9) 


/ 9 

IT 1 r 

—n 


P/ 2 2 -2 ^ 9 - —9 


m 


“ 2 a. {(r -x„ i5c ^4 . ) 
1 n-i-l n+l-i' 


2 - -9 

m 


where g(x) is as in (8). 

2 A 

Numerically computed results for L ( c^) -optimal 

S-stop methods are presented in Tables 5,1 and 5, 1(a) -(d), 

Toblc 5,1 displays the coefficients and the error norms 

for the optimal method. The optimal error norms at the 

beginning of the table have now reduced further (as 

canpared to Chapters 2-4) , The general behaviour of the 

coefficients however# has been observed to be as in the 

earlier diapters. As seen from Tables 5. 1(a) -(d) the 

oxislicit L^(c ) -optimal multistep method has been just better 
XT 
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on equations 13/ 16 and 24; one decimal place better on 

equations 2/8/10/12,15/17,18/19, 20 / 21 , 22 and 23^ two 

decimal places better on equations 1,3, 4, 5, 6 and li; and 

three decimal places better on equations 7 and 9 / in 

comparison to the 5-step Adams -Bashforth method. As 

compared to 5-step Adams-Moulton method, the implicit 

,/s 

L 2 ( 0 ^) -optimal 5-step method has been just better on 
equations 16,15,17,18 and 24; one decimal place better 
on equations 5,10,11,12,19/20,21,22 and 23j two decimal 
places better on equations 1,2,4 and 7>and three decimal 
places better on equations 3 , 6,8 and 9, 


5.4 


Optimal Multi step Methods in L^(c ) Interpol atory 
for Polyncmials 


2 

As was done in Section 3.6 for the space H (c^), we 
consider here, for the space L^(c^), the optimal multi step 
form\ilae interpolatory for polynomials corresponding to 
the special situation of equispatsd points on a line within 
G^, Let the distance between the consecutive points be h. 

■^P 2 ^ s 

The coefficients 21 of an L (c^) -optimal mxxltistep 
formula 


( 10 ) 


m 

yn.l- 


p -1 

s 

j =0 

m-6 


’tl'^n+6 


tl 


to 


J=P 




■"tl “’^^tl 


interpolatory for polynomials of degree p<m't 6 ^^/with 



156 


prefixed a.‘s such that Sa. =1 and with r. ’s as for 
the corresponding usual formiilae, are to satisfy the 
normal equations 


( 11 ) 

where 

A p 

(12) = 


c r 


d^ 


-^p /'p 

* * "im-6 


T 

) . 


to 


(13) c?.= 

' —ij 


71 


K=0 1.0 (r 


(14) ^ 




-3 


m 


q=l ^ 




2r2 P-1 „ q 1 H.V i a 

■— s y^ { 2 2 (-1)^ ^(^) (S ‘ 

^ q=0 h=0 1=0 ^ ^ 


r i- 


^+ 6 ^ 1 - 3 < ^ n + 6^^-1 


} , 


i-fj = p(l)m—6 


to 


The above described method is characterized as 
follows : 

Theorem 3t The L(c^) -optimal multi step formula (10) 
interpolatory for polynomials of degree p < is 
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characterized by being interpolatory for functions 


h. (x) = S 

J k=0 ^ 


X 




, j=^(l)m-6 


to* 


Proof ; The proof is an easy consequence of (2) and (11), 

y^'p 

TO compute bj, j = 5^Q(l)m/ in the alternate 
representation 

< yntl 


m 


m 


i=l 


of (10)/ we use the relations 


h .S bj yn+l-j> 


J=6 


to 


(16) 

—l+o . _ 
to 

^ m-6 
= (-1)-^ S 

jsSl 

where 



(17) 

II 

il? ' p . 


.i; 


to 


"to 


It is important to note that ^ different in the 

ejiplicit and the implicit situations, 

P 

In Tables 5, 2-5, 6 we present the coefficients bj 

2 

and the error norms of the L (c^) -optimal 5-step methods 
interpolatory for polynomials of degree p = 1(1)5. We 
observe that the general behaviour of the error norms and 
the coefficients is the same as in the case of Chapters 2-4. 

A summary of numerical results obtained by using these 
coefficients on the set of twenty four differential equations 


is as follows; 
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As compared to a 5~step Adams-Bashforth method^, the 
ejqolxcit L ( c^) -optimal method with p=l has performed 
just better on equations 17/18 and 24? one decimal place 
better on equations 2/ 6/3/10/11/12/19 / 20 , 22 and 23, tv/o 
decimal places better on equations 1/3/4/5/9 and 21;and 
three decimal places better on equation 7, On equations 
13/ 14' and 16/ however/ the optimal explicit method has 
been just worse than the usual rf.athod and conparable 
on equation 15. In comparison to 5 -step Adams Moulton 
method the implicit L (c^^.) -optimal method with p=l has 
been comparable on equations 14,15 and 16; just better on 
equations 13/17/18 and 23/ one decimal place better on 
equations 5/10/12/19,20/ 21/22 and 24x and two decimal places 
better on equations 1/ 2, 3/4/ 6, 7/8/9 and 11. 

The optimal explicit method with p=2 has performed 
just better on equations 17/18 and 24; one decimal place 
better on equations 2/8,15,16/19/ 20/21/ 22 and 23/ two 
decimal places better on equations 1,3/ 4, 5, 6,9 and 10/ and 
three decimal places better on equation 7. The implicit 
method with p=2 has been just better on equations 15,16,17/ 
18 and 24, one decimal place better on equations 5,10/19/20/ 
21/ 22 and 23/ and two decimal places better on equations 
1/2/ 3/4/ 6/7,8 and 9. Both the optimal methods have hsen j 
comparable on the remaining equations 11,12/13 and 14, j 
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The optimal explicit method with p=3 has performed 
just better on equations 15 and 24J one decimal place 
better on equations 1/2,3/5,6,8^9,10,19^22 and 23; and two 
decimal places better on equations 4,7,20 and 21, The 
optimal implicit method has been just better on equations 
15 and 24; one decimal place better on equations 1,2, 4,5, 
10,19, 20, 21, 22 and 23 ! two decimal places better on 
equations 3,7,8 and 9; and three decimal places better on 
equation 6, The implicit and explicit optimal methods, 
both, have been comparable on equations 11 and 16* and just 
worse on equations 12,13,14,17 and 18 as compared to the 
usual methods. 

The optimal explicit method with p=4 has performed 
just better on equations 2,10 and 24? one decimal place 
better on equations 1,3,4,5,6,8,9,19,20,21,22 and 23 ? and 
two decimal places better on equation 7. The optimal 
implicit method has been just better on equations 10 and 
24? one decimal place better on equations 1, 2, 3,4, 5, 6,7, 
8,19,20,21,22 and 23j and two decimal places better on 
equation 9, Both the implicit and ejpilicit optimal methods 
have been comparable on equations 13,14,15,16 and 17j and 
just worse on equations 11,12 and 18, 

The implicit optimal method with p=5 was found just 
better on equations 10,15,17; one decimal place better on 
equations 1,2,4,5,8,19,20,21,22,23 and 24? and two decimal 
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places better on equations and 9. This method 

has been comparable on equations 13/ 14 and 16; and just 
worse on equations 11,12 and 18. 


2 ^ 

5,5 Optimal Multi step Methods in L (c^) Interpol atory for 
a Set of Preassiqned Functions 

The coefficients bj of an optimal multistep forraula 
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in L^(c^), interpolatory for functions { 
con be obtained by solving linear equations 
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where "c, d are as in (4), X = [ * 
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and 0 is a p x p noill matrix. 

The above method is characterized by 

Theorem! 4: An’ optimal multi step method in L^(c^) inter- 
polatory for linearly independent functions ^ i 

is locally interpolatoiry for functions 
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X 


m 
2 


X 


(r2-x„^^..x) 


k = 5^Q(l)m-p , 


where 


{^%q^i=l(l)p = p"^E, 

q=l(l)m-p+6^^ ^ 

in which P = [ 9.' (x_^^ .)]. ./-v is assimed non- 

'- 1 n'+'l-j -‘i=l(l;p 

j=m+l~p(l)ra 

singular and E = [ ?! (=<nH.i.j)]i=l( l)p 

Proof I The proof is immediate from Theorem 1,8, using (2), 

In Table 5,7 we present the error norms and coeffi*- 
2 

cients of L (c^) -optimal 5-step methods interpolatory for 
fmctions 9^(x) = exp(l,6x) and 92(x) = exp(-1.6x). The 
genercil behaviour of the coefficients and the error norms 
in the present case is found to be the same as is the 
case of Chapteis 2 to 4, The performance of these methods 
as compared to usual 5-step methods on the twenty four 
differential equations has been as follows; 

The explicit optimal method interpolatory for 
exp(1.6x) and exp(— 1.6x) has performed just better on 
equations 16,17,18 and 24; one decimal place better on 
equations 15,19,20,21,22 and 23; two decimal places better 
on equations 1,2,3,4,5,6,8,9,10; and three decimal places 
better on equation 7. In the iit^licit case the method 
has been just better on equations 17,18,23 and 24; one 


8 
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decimal place better on equations 5,15^16,19/20^21 and 22 
two decimal places better on equations 2,3, 6,7,8,9,10? 
and three decimal places better on equations 1 and 4. 

The e3q>licit and implicit methods^ both, have worked compa- 
rably on the remaining equations 11 , 12,13 and 14. 


5.6 Limiting Behaviour o£ coe f ficients as r «» <» 

In this section we describe the limiting behaviour 
of the coefficients of optimal end various interpolatory 
optimal multistep formixLae in the space L (c^) • Since 
the remarks made at the beginning of Section 4*6, ccaaceming 
the proofs of various results, remain applicable in the 
present case as well, only the statements of the restolts 
are given. 

Theorem S: Let 


( 20 ) 


^n+l 


m 

2 

is:l 


^i^n-i-'l- 


m 
h S 
j=^' 




to 


be a multistep formTola interpolatory for functions 
’’m’ and let the matrix [q 1, 

be non —singular. Let 


m m ^ „ 

1-:*1 1 = 6 , 

j ^-to 


2,/> 


be an optimal multistep formula in L (c ) interpolatory 


6 ^ ° ^ “"^to 

to ^ ^ to 


for functions { f 


•tki 
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(p=0 being the non-interpolatory case). If the local 
truncation error fxonctional for (20) satisfies 


./\P 

=0 / J = p+6.^(l)m , 


JC *^oo 


—nr j 


to 


then 


lim = b^ 

x: oo J 


/ j = 6^Q(l)m. 


Theorem 6, Let a^, i=l(l)m be given constants satis- 

2 ^ 

fying Sa^= 1. Then the L (c^) -coefficients bj^ of an 
optimal multi step formula 


( 22 ) 


^n+l 


isil 


“i^n+l-r 




approach the coefficients of the corresponding usual 
formula as r oo 

Theorem 7i Let a^^# i=l(l)m,be given constants such 

2 A p 

that Sa. = 1. Then the L (c ) -coefficients b. of 
i r 

an optimal milLtistep formula 


(23) 


Y. 


n+1 


m 

S 

i^l 


^i^nfl-i"** 


h 


m 

E 

j=5 


A p 
— jr 


'^n+l-i^' 


to 


interpolatory for polynomials of degree p < m "*■ ^tl' 
approach the coefficients of the corresponding usual 
formula as r -* <» . 

Theorem 8; Let a^^, i=l(l)m, be given constants such 

that Sa. = 1. 

1 


Let the matrix 
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M = 


^ ( X ) 




n+6 


tl 


m-p~6 


X. 


to 


■n+6 


tl 




’4Kfl-ra> 


X 


'nihl-ina 


X 




'n+l-m 


2/A . . . A p 

be non -singular. Then the L (c^) -coefficients bj^ 
of an optimal multi step formula 


m 


m 


(24) 


1-1 -J'^to 


interpolatory for the functions { p < 

P 

approach the coefficients bj of the unique multi step 
formula 


m 


m 


(25) 


2 m-pt6^2. 

interpolatory for functions .../9 , x,x 5. 


5,7 A Comparison of Different Methods 

In Table 5,8, the end point (x^^=1.2) errors for 
the set of twenty four differential ecguations for 1 to 4 
step ej<plicit and implicit L^(c^) -optimal methods , along 



Table 5*8 


situation; y in l cc r) t Rc 2.0l H* #10 

ERRORS AT END POINT FOR EXPLICIT METHODS 


Eq.NO, 


USUAL OPTIMAL USUAL OPTIMAL 


OPTIMAL 


Optimal 


-0,41E-0l 

-0,32E-0J 

-D.50E-01 

-0.39E-01 

O.UE^OO 

-0.44E^52 

-0.50E-03 

-0,81E-62 

-0.20E-02 

-0.17E+01 

-0.13E+54 

0.13E-f04 

-D,99E-02 

--0.85E-02 

-0,80E-02 

-0.69E-02 

-0^67E-05 

-0.66E-04 

0.63E-01 

-0.60E*02 

-0*57E-02 

0.22E-01 

0.26E+00 

0.97E-04 


-0,18E-01 

-0.13E-01 

-0,25E-5i 

-0.19E-01 

0,44E-01 

-0,24E-02 

-0,3lE-03 

-0.42E-02 

•0.12E-02 

-0.11E401 

-0.36E403 

0l34E’^03 

-0.30E-02 

-0.2iE-02 

•0.27E-02 

-0,23E-02 

-0,39E-05 

-0.38E-04 

-0,21E*01 

-a.i2E-oi 

-0j3E-01 

-oa6E-ai 

0.47E-01 

-0.59E-02 


0*25E-02 

8;iS|:S 

0.28E-02 

-0,46e-02 

0.57E-03 

0.96E-04 

olltiloi 

8;Ue+03 

-0.77E+03 

0.18E-03 

0,14E-03 

0.20E-03 

0.16E-03 

0.82E-06 

0.79E-05 

0.53E-02 

0.54E-03 

0.72E-03 

0.28E-02 

-0.36E-02 

-0,65E-04 


8:iii:8i : 

-P.62E-03 
-0.30E-04 - 


-0,13E-Ol 

-0.46E+Q2 

0,37e+02 

0.73E-04 

0.71E-04 

0.30E-04 

0.32E-04 

-0.24E-06 

-0.24E-0S 

-0.52E-02 

-0.84E-03 

-0.11E“02 

-O. 36 E -02 

-0.75E-02 

-0.52E-03 


0.30E-02 

8:Ui:8i 

8:181:81 

0.69E-03 

O.llE-03 

0.94E-03 

0,37E-03 

0.21E+00 

0.29E+03 

0.29E+03 

0.l7e-03 

8:18i:8i 

8:111:81 

0,69E-05 

0.71E-02 

knr-u 

8:111:81 

a.77E-04 


•0.65E-03 

0a6E-02 

•0.80E-04 

•0.93E-05 

:8:ir-8! 

• 0 . 36 E -61 

•0.98E+51 

0.62E+01 

-0.12E-03 

•O.lOE-03 

•0.97E-04 

■0.83E-04 

•0,26E-07 

■0.24E-06 

0.59E-03 

-Q.27E-03 

•0.16E-03 

8:iS|:8l 

'»0*12E»03 


8:?1i:81 

■klilzU 

8:l^^:81 

MilzU 

B »E^U2 
8E405 
0*23E-04 

0*3Yf04 

7JiE-»04 

e Ofc.-06 
iL-06 
0.4DE-a2 
a,7/E-o5 

*Ol20E-02 

Ols/E-Oi 


-0.i6E-04 

•U.20E-04 

-0.39E-55 

-5,22E-54 

0,65e-04 

-U,59E-05 

a.l6E-05 

•0,13E-04 

0,47E-05 

-0:i5E-02 

-0*10E+01 

-0,69E+00 

-0,49E-05 

-g*41E-05 

•U*46E-05 

-U,39E-05 

0,18E-07 

-8:i?|:8l 

:8:f6i:88 

-0,75E-03 

-0,UE-02 

•0,22E-03 


ERRORS AT END POINT FOR IMPLICIT METHODS 


Eg. No, 


USUAL optimal 


OPTIMAL USUAL OPTIMAL USUAL OPTIMAL 


-0*19E-03 

-0,l2E-03 

•0*31E-03 

-0.20E-03 

0.3eE-03 

-a.45E-04 

-0.79E-05 

-0*68E-04 

-0.28E-04 

-O.llE-01 

-0.18E'»'03 

0.18E403 

-0.15E-04 

-0a2E-04 

-0*16E-04 

-0.13E-04 

-a.73E-07 

-0,71E-06 

-0,28E*03 

-0.63E-04 

-0.79E-04 

-0.17E-03 

0.32E-03 

0.46E-0S 


•U*31E*04 

- 0 I 18 E-O 4 

-0,50E-04 

-gl28E-04 

5*87E*04 

-0.25E-05 

-0.52E-06 

0,22E-05 

-0.29E-06 

-0*40E-03 

0.53E+0J 

-0.25E401 

•0.66E-a5 

-0*62E-05 

-0.36E-05 

-0,34E-05 

0.14E-07 

0.14E-06 

0,42E-03 

0.94E-04 

0.12E-03 

0*25E-03 

O.71E-03 

0 I 66 E-O 4 


0*27e-03 

0.16E-03 

0,58E-04 

0.94E-05 

0.84E-04 

0.33E-04 

0.21E-01 

0.28E+02 

-0.28E+^02 

O.ieE-04 

0.14E-04 

0.20E-04 

-O.ilE-04 


0.78E-04 
0.55E-04 
0.12E-03 
0.78E-54 
-0.20E-M 
0.89E-05 
0.88E-06 
0.12E-04 
0.22E-05 
0.41E-02 
0.14E-0i 
0.S3E'i'OO 
0.J5E-04 
0.13E-04 
0.12E-04 
O.lOE-04 
-0.16E-08 
-0.19E-07 
-0.47E-04 
0.41E>04 
0.23E-04 
-0,36E-04 
-0.49E-0| 
0.19E-04 - 


-D.58E-04 

•5,28E-04 

-5.12E-03 

-0.526-04 

O.eOE-04 

-0.17E-04 

-0.28E-05 

-0.19E-04 

-0.76E-05 

-0.25E-02 

-0.60E+01 

O.eOE't^Ol 

-O.lOE-05 

-0.68E-06 

-0.14E-05 

-0:96E-06 

-0.46E-08 

-0.36E-07 

-0.13E-03 

-0.25E-04 

-0.23E-04 

-0.55E-04 

0.72E-04 

-0.19E-05 


0,92E-0§ 

0.13E-05 

•O.3OE-06 

.8:886:81 

O.llE-07 

-0.14E-06 

.8:tll:8t 

0.64E-04 

■8:881188 

0.54E-06 

0.49E-06 

0.41E-06 

0.36E-06 

-0.17E-08 

-O.ieE-07 

0.88E-04 

0.20E-04 

0.19E-04 

0.39E-04 

0.49E-04 

0.12E-04 


0,6SE-04 

0,2»E-04 

0.14E-03 

o:6t>e-04 

-olsze-ol 

0.21E-04 

0,34E-05 

O.24E-04 

0:8bE-05 

0 . 2 oe -02 

O.lbE+01 

-O.lbE+Ol 

0.69E-06 

0.43E-06 

O. 9 /E-O 6 

0.6BE-06 

0:3UE-08 

(),2&E-07 

6:ibE-03 

0:30E-04 

0:2SE-04 

0.60E-04 

-o:7*E-04 

6:i»e-05 


0.44E-05 

0.30E-05 

0.|9|-05 
-0.126-04 
0:481-06 
0.57E-07 
0. 486-06 
0:iOE-06 
0,15E-03 
-0.59E-01 
0.13E+00 
O.lOE-05 
0.e9E-06 
0.74E-06 
0.65E-06 
-0.25E-09 
-0.14E-08 
-0.27E-04 
-0.|7E-05 
-0.376-05 
-0.13E-04 
-O.47E-04 
-0.18E-05 
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with those for the usual methods, are presented. We / 
observe that the optimal methods have performed' better/ 

I 

in all the cases except in the case of equation 24. \ 

A 

In Figure 5.1 we have plotted log^^C t IT^ 1 1 ) and 
logic ( I 1 I ) different points for explicit methods 
of 1 to 5 steps. In Figure 5,2 the same is done for the 
implicit methods. The curves in these figures indicate 
the same general behaviour of the local truncation error 
norms as in the earlier chapters. 

In Table 5,9 we present the end point (x^=1.2) 
errors in the numerical solution of the twenty four 

/V 

differential equations by five step usual, optimal (M), 

jL 2 

optimal interpolatory for polynomials vM. ' /.,.#etc., 

the superscript denoting the degree of the polynomials) 
and optimal interpolatory for exp (i.6x), exp(-1.6x) 

/\ P 

(denoted as M )-methods, in the es^licit case. We 
observe that for these explicit methods the error is 
least for the method: "m on equations 11,12,17,18,19,21 
and 22j M"*" on equation 9^ on 1,3/4 and 6^ on 

A A 'A.P , 

equations 20 and 23? M on equation 24} and M on equations 
2,5,7,8,10,13,14,15 and 16. 

In Table 5,10 we, similarly, present the end point 

A A.1 A .5 

2) errors for the implicit methods and 

We observe that the error is least for: M on equations 
3/ 6, 7,8, and 9; on 11,12,19,21 and 22} on equations 
18,20,24; on equation 17; on equationsl3 and 23; 





Figure 5.2: Local Tnincation Error in Implicit L (Cj.) 

Optimal Multistep Methods/Usual Methods 

Cr=2.01/ h=0.l) 




Table 5.9 

End point 2 errors for e3<plicit L^(c ) -optimal 

Multi step "nethods ^ 
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CHAPTER 6 


ERROR ANALYSIS OF OPTIMAL MULTI STEP METHODS 

6,1 I ntroducti on 

This chapter is a study of ceri:ain theoretical 
aspects of the optimal mxiltistep methods discussed in the 
previous chapters. These aspects concern convergence/ 
magnitude of discretization error and propagation of 
round offs in the implementation of the optimal and 
quadrature optimal mxiltistep methods. The study also 
includes the optimal and quadrature optimal methods which 
are interpolatory for polynomials of a certain degree. 

The results obtained in this chapter establish the 
applicability of our methods to a general situation in 
which the solution does not necessarily belong to the 
underlying Hilbert space. 

In Section 6,2 we prove that the coefficients of 
optimal mxoltistep methods as well as those of quadratxire 
optimal multistep methods approach the coefficients of 

the^- related usual 'methodd 'as h 0, ' al'So if th^ 

optimal/quadrature optimal methods are interpolatory for 
polynomials of certain degree. In section 6,3 we state 
sane results fron Henrici [27] which are required in the , 

subsequent sections. In Section 6.4 we give an estimate ‘ 
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for local triincation error. In Section 6.5 convergence 
of these methods is established/ while in Section 6,6 
a bound on discretization error and in section 6,7 a 
bound on round off error propagation have been derived. 

Throughout this chapter, H denotes either of the 
spaces H Cep or b (cp as described in chapters 2 
and 4. The points Xj^*s satisfying a=x^<x^<,, ,<3<^j=b 
are assumed to be equispaced with spacing h. The interval 
[a/b] of integration is a subset of the open interval 
(-r/r). An m-step method would then start from the point 
and will calculate the solution at the points x^ 
through x^# Note that N-is a fuinction of h so that 
Nh = b— a. 

In the sequel an H-optimal multi step method 
commonly denotes any of optimal/ quadrature optimal/ 
interpolatory for polynomials of specified degree methods, 
implicit or explicit. 

Our main concern in the earlier chapters has been 
the computation of coefficients of multi step formulae 
for a fixed set of nodes. In this chapter, however, we 
need to stress the dependence of the coefficients bj's 
on the current variable point integration 

scheme. For this reason the dependence of bps on the 
point is indicated by writting them as 

In section 6.3 onwards, where a slight change of notations 


I 
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occurs/ the dependence of coefficients /3j's on the point 
has been denoted by presence of a siabscript m in 
where 5?^ is the farthest point from the current point 
in a k step formiola. 


6.2 Macmitude of the Coefficients in H-optimal 
Multi step Methods 

The following two lonmas relate the coefficients 
in a H~optimal raultistep formula with the coefficients 
of the corresponding usual formula. 


Lemma 1. Let 
( 1 ) 


m 

to 


be an m— step quadrat-ure optimal multi step methocV^inter- 
polatory for polynomials of degree p < m+ in H and 
let 


m 


( 2 ) 




JS:6, 


to 


denote the corresponding visual mviltistep method. Then/ 
uniformly for m-l^nj^-l. 


'^j/n+l“^j' = 0(h)/ h - 0/ j = 6^o(l)m. 


Proof : The methods (1) and (2) are implicit if t=l 

and explicit if t=0. Correspondingly let us consider 
polynomials 
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(3) 




fr “ ^^ 1 -i 

jA 


ni-l<n<N-l^ 5 <l<m 

»— • 13.0“ “ 


As ® belong to [a^b] and they are equi spaced, the 

norm of the polynomial L. in H satisfies 


(4) 




-(m-6 

0(h 



/ 


uniformly in i and n. Since the formula (2) is exact 

for polynanials of degree m’^*l“6^Q we have 

x: 1 

hbj » / Lj^(x)dx, j«6^Q(l)m. 

^-s 

Jipplying the method (1) on Lij^(x) and using lemma 5 
of the next section (Lemma 5.7/ Henri ci [27]) 




I / 


x„ 


Ljj^(x)dx 


‘inti’ i "V" 


On" 


n^s 


< iiLjj,! 1 n(^ II 

m'«-2-6 

< C. LIL. II h 
— t jn 


. . ,s 
^‘*‘^"^to 


= 0{h^) 


uniformly in n where/ in view of Theorems 2*1 and 4.1/ 

rati- 6 


^rr,^T = ’^ax I ID^ 

X e [a,b] ^ 


•^to 


1 1 < 


oo / 


m+1-6, 

D denotes the derivative evaluation functional of 

X 

order irH-1-6^^ and ^ and designate the quadrature 
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error functionals in the optimal and usual cases/ 
respectively, corresponding to the point Hence 

(5) 1 = 0(h) / h - 0 , 

uniformly in j = 6^Q(l)m/ m~l < n < N-l, completing the 
proof. 

Lemma 2: Let 

m m ^ 

i^=:l 

m-l < n < N-1 

be an optimal mxiltistep raetho<^interpolatory for polyncmials 
of degree p < m 6^^ in H and let 

ra m 

1*1 J* to 

m-1 < r ;C N— 1 f 

be the corresponding usual method. Then, uniformly for 
m~l < n < N-1, 

‘-^j/h+l “ h ^ 0, j = S^^CDm. 

Proof : As before t = 1 refers to an inplicit and t = 0 
to an e^^plicdt case. Let 

X 

(8) ^ ^to ^ ^ ^ S ^ S 

a 

where L. (s) is as defined in (3). It follows from (4) 


that 
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(9) ^ 

uniformly in i and n, since the usual formula (7) is 
exact for polynomials of degree m+6^^/ applying it on 
Min(x) we get 

m 


Next# by applying the optimal formula (5) 
and reasoning as in the proof of Lemma 1# 

m 





S 

i^il 


^i^-n^^n^l-i^ 


on ix) 

we get 


/V 

- h b, 


J/n+] 


< IIM-J^ll lITJI 

< I IM.^I I I I T^l I 

i c° say 

= oCh^) 

Toniformly in m-1 < n £ N-1. Here 


m'''l-6. 


M 


m+1-6 


max 


to X e [a#b] 


HD, 


to 


X 


1 1 


[a/b] is the interval of integration and and 

denote the local truncation error fxonctionals corresponding 
to the point optimal and usual cases# respectively* 

Thus# we have 


( 10 ) 




uniformly in m-1 < n < N-1 


OCh)# 
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6.3 Some Auxlliarv Results 

In this section we state three results from 
Henrici [27] which are to be used in our later analysis. 

Let 

(11) L(y(x);h) = (Xj^y(x+kh) + 

(x+}ch)+j3j^_^y* (x+(k-l)h)+. . ( x) } 

be a difference operator. The operator L(y(x);h) 
corresponds to a general k- step method given by 

where 

W = w> ' ^ • 

Here^ if 0 the method is explicit and it is inplicit 

otherwise. 

The difference operator L, or the method (12)^ is 
said to be of order p if 
k 

(13) 2 a. = 0, 

i;a0 ^ 

k 

(14) 2 i^ a. = 

iisl ^ 

and 

(15) E a. 

1^1 ^ 


q s J"* ^ . q = l(l)p» 

J.l ^ 


jifpHlS jP Pj . 
J -1 ^ 
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A raultistep method given by (12) is called 
consistent if 


(16) 

p(l) 

= 0# 

and 

(17) 

P' (1) 

= 0(1)# 

where 

(18) 

P(0 


and 

(19) 

o(i) 



are the associated polynomials of miiltistep formvila (12). 
The conditions (16) —(IT) ernount to saying that the formula 
(12) is locally exact for pol 3 mamials of degree one# or 
that the order p of (12) is at least one. 

The formula (12)# or the method associated vd.th 
it# Tfthere prefixed# is 'usual' if the 

B ' s are choosen so as to maximize the order of the 
method. Note that for a xisual k-step method the order 
is at least k in the explicit case and k+1 in the 
implicit case. 

1 

Lemma 3J Let the polynomial p(0 = * ^-1^ +...+% 

satisfy the condition of stability and let the coefficients 
y^(l=0#l#2# .. .) be defined by 

i — s 
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Then, 

(20) r = sup ir 1 < CO . 

1 = 0 , 1 /... ^ 

Lemma 4. Let the polynomial P(s) = +. ..+a^ satisfy 

the condition of stahility, let B"''', /3 and A be nonnegative 
constante such that 


(21) ‘^0,n^ ^ 

I S. ^ f ^ ^ ^ ^ f n=sO/ 1/ 2, • • . ,N 


■.-1 


and let 0 < h < . Then, every solution of 


(22> “kW ♦ Vm 


for which 



(23) 



satisfies 



(24) 

\z 1 
n 

, nhli « ^ 

^ ^ 4r 

where 



(25) 

L* 

= r*B* , K* = r* (N A + AZk) , 

(26) 

A 

= JcXo‘ 

and 

r 

r 

i+hioj^r^^ * 

Lemma 5. 

Let 

L (y(x)jh) be a difference operator of 

order p>0 

• There exists a constant G>0 depending 
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only on L/ such that 

(27) |L(y(x) ; h) I < h^**"^ GY acx/ xt]<h^ 

for all fianctions y(x) having a continuous derivative 

of order in [a/b] with Y = max J (x)I . 

xe [a/b] 

6.4 A Bovind for Local Truncation Error 

A generalization of Lemma 5/ applicable to 
H-optimal methods, is as follows: 

Lemma 6, Let L(y(x)jh) given by (11) be the difference 
operator of a 3c-step usual method of order p and let 

(28) L°(y(x)jh) = a^y(x^kh) + a^_^y(x+(]<;“l)h) . .+aQy(x) 

- h{i3^^^' {x+kh)t, . (x) } 

denote the difference operator of a corresponding 
H-optimal k-step method (the subscript x in ^ denoting 
the dependence of these coefficients on x) . Let y possess 
a continuous derivative of order p+l on the interval 
[a,b] and set 

(29) Y= max ly^^^ (x) 1 

X e [a/b 1 
0<i<p*fl 


Then 
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(30) lOj^yCx+kh) aj^_^y(x+(k-l)h)+. , .•♦■a^yCx) 

“ ^^^o,x Y‘ y'(x)} I 
= lL°(y(x);h)l < hP'*‘^G*Y^ 

uniformly in X/ x+kh e [a/t)"] where G'‘ is a constant 
independent of y. 

Proof ; We know that 

yfh ^ 

C 3 I) y(t) ■= y(x) + (t-x)y' (x)+...» r J- (t-s)Pyfey' ds 

JP # jP » 

= P (t) Q(t) , say. 

Then# if (26) is an optimal method/interpolatory for 
polynomials of a certain degree (as distinguished from a 
quadrature optimal case) and liL'^Il# 1 IL II respectively 
denote the norms of the functionals and L corresponding 

to the current point x-ekh, we have 


(32) lL^(P(x)jh)l <IIPII IlL^Il 

X X 

< 1 IPI 1 1 IL I 1 

“ X 

< (C3_y) (C2Mp^3_hP'*‘^) # 

using Lemma 5 and (29) (for estimating 1 (Pit)# where 



max 

X £ [a#b] ^ 




and 


are constants not depending on x 
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By directly substituting the expression for Q(x) in ( 28) 
and simplifing we get 

(33) lL°(Q(x))h)l < C^Y , 

where is a constant independent of x. Prom (31)/ 

(32) and ( 33 ) it follows that 

(34) lL°(y(x)jh)l < h^"^^ O'Y./ 
where we have set 

G* . (0,C2» C3)Mp^^. 

This proves the result in the optimal case* 

If the method corresponding to the operator 

in (28) represents a quadratiire optimal method/interpo- 

latory for polynomials of a certain degree/ rewriting 
o ^ 

L° and L tenporarily as and L*?, respectively/ 

•X. X. (X X. 

and regarding these as quadrature error functionals 
corresponding to the current point x+kh/ we have 

(35) lL°(P(x)»h)l <1IP‘11 

-X yC 

< UP' 1 1 1 ilJi I 

< (C^Y) (C 5 hP’^^p) 

where and are constants independent of X/ and 

M = max I IdJi I . 

P xe [a/b] 
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AlsO/ by a direct sxibstitution, we have 

(36) lL°(Q(x);h)I <C.YhP"'^ , 

the constant Cg being independent of x* Combining (3l), 
( 35 ) cjnd ( 36 )^ we get 

(37) lL°(y(x)5h)l < h^* V y , 
where we have put G* = (C.C + Cg)M , 

5 o p 

Hence, from (34) and (37), the result is valid 
for any H-optimal method, 

6,5 Convergence of H-optimal Methods 

Let f(x,y) satisfy the following conditions: 

(a) f(x,y) is defined and continuous in the strip 
a-<iX<b, - °°<y<»®, where a,b are finite, 

(B) For any two numbers y and y^ and xe [a,b]J ^ a 
constant L exists -stich that 

If(x,y) - f(x,y*)l <L!y-^I. 

Then, a linear m\iltistep method 

‘38) Vn-ric’- Vl^n+k-1-^”--" Vn 

where f^ = f(Xjn*^m^ (m^O, 1, 2, . . , ) , with ^ 0 and 
la^l +• > 0, is called convergent if for all 
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fimctions f(x#y) satisfying the condLtions (A) and (B) 
and for all values of r? there holds 

(39) lim y = y(x) 

h >^0 ^ 

x„=x 

n 

for all xe [a/b] and all solutions £7";^^ diff- 

erence equation (38) having starting values y = '^L,(b) 
satisfying 

(40) lim r) (h) = 7) / M = 

h **0 ^ 

where y(x) denotes the solution of the initial value 
problem 

(41) ^ = f(x,y) , y(a) = '0 . 

The criterion of convergence described above is 
similar to that given by Henrici [27, p. 217]. 

We now prove. 

Theorem 1: An H— optimal method corresponding to a stable 
and consistent usual m^ultistep method is convergent. 

Proof . Let 

(42) 6 = 5(h) = max II "H (h) -y( a+jah) I 

with 

lim 6(h) - 0. 

h-*0 

Let the difference operators in the usual and H-optimal 
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H 

cases be denoted by L(y(x);h) and b^(y(x)5h)/ 

respectively. In view of Lemmas 1-2 we have 

I 1 = 0 Ch) where coefficients 

of the H-optimal method occurring in L^(y(x);h) 

corresponding to x = With i ly'l I = nmx ly* (x) 1, 

X e [ 

we have 


(43) 


lL^(y(x)jh) -L(y(x);h)l 

jsaO 

= Oih^), 


JlyMi 


unifomly in x e [a/b-kh]. since L is consistent and 
y* = f(x/y(x)) is continuous in [a/b]/ 


(44) lL(y(x);h) 1 =o(h), 
uniformly in xe I^a#b-kh]. Therefore 

(45) iL^(y(x);h)I < he(h) 


where 9(h) ** 0 as h ^ 0 and is independent of 
X e [a/b-kh] . Thus with 


(46) 
and 

(47) 


®m “ ^ ' m=0,l#2, ... 




'’ra 


0 / if e = 0, 
m 


we have 
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(48) 


O; '3 . .+ 0 C e 

Jc m-fk o m 


= e^h ©Ch) 

where 1©^I < 1, The defini-tion (47) of g_/ in view of 
m — 

the Lipschitz condition (B) as stated in the definition 
of convergence/ implies that 


^ m - 0/1/2/... 

/^plying Lemma 4/ with z^= e^^z Z = 5 (h) / A = h© (h) / 

N = (x — a)/h and B*s= BL where 
n 


(49) B . sup (If^l . lO*— ^ '4n'5 < “ ' 
m/h 

(by Letimas 1- 2)/ we have 

< r [(x^-a) ©(h) +A6(h)k] e3?p[ (Xj^-a) P B*] 

where 

p* _ P 

l-rtilo^P^ /3L 

with = sup Lonmas 1-2)/ A 

as in (26) and P as in (20). Since ©(h) 0 and 

6 (h) •♦0 as h 0 we have 1 e^^I ** 0 
completes the proof. 


as h 0. This 
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^•6 A Boimd for the Discretization Error 

Let us assxame that the exact solution y(x) has 
a continuous derivative of order p**"l on [a/b]/ with 

P+1 n-) 

(50) Y= max S Jy''^''(x)l. 

xe [a/b]is0 

With the milder assxamption 

< 5 ^ > “kW “k-l W-l*- • • Vn 

where and q are non -negative constants and <, 1 

(rather than the exact satisfaction of the difference 
scheme) / we have the following a priori boxind on the 
discretization errors 

^ —1 

Theorem 2t If h is such that hl^j^ Oj^ 1 L < 1/ where 
jS* = max and x^e [a^b], 

^ r*[A6h + (Xj^-a) (K^h^l+G^Yh^)] exp (x^-a)L b] 

where 

( 52 ) r* 

l-blOj^^ /^IL 

and 6 t A,G^f P and B are as in (42), (26), (30), (20) 
and (49), respectively. 

Proof . Subtracting L^ fecm (51) we have 
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= g"'^ hP’^^l , m=0,l,2,.., 

by Lemma 6, where < !• 1'he theoron now follows by 

an application of Lemma 4. 

Note that if 6 (h) = O(h^) and q > p then e^= O(h^) . 

6, 7 An a Priori Bound on Round off Errors 

To study the influence of local round off errors, 
instead of the exact difference scheme 

Vn^k'^- • •♦“c.yn" ^ Vk'^n-rt^’ ‘ ' 

<v> 

let the quantities y^^, actually calculated, satisfy 

^n+k' 2,... 

where the quantities are the local round off errors 
with 

l^n+k' ^ ^ / n = 0,1,2,*.., 

where e is a constant. 

An a priori bound for the accumiiLated round off 

errors r = y -y is, then, given by the following: 
n -'n -^n 

Theorem 3i With h and ^ as in Theorem 2, 

•-*'1 «». " 5 ^ t “ T 

^ ^ * expLCx^^a) 1 BLJ , n=*, k+1, ... 
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where i is as in (52)# B is as in (49) and L is the 

\ 

Lipschite constant for the function f(x#y) in (B). 

Proof . Along the lines of the proof of Theorem 2, if 

9m “i 

L 0 . If 0, 

we have 

The result follows from this by an application of Lenma 4# 
with z^= r^^/ A = e /N=(:!^-a)/h and Z=0 (i.e. r^=r^=. . .=r^__^=0>. 

It thus follows from Theorems 1-3 that the H-optimal 
multi st^ methods of the thesis continue to enjoy the 
approximation properties of the corresponding consistent 
and stable usual methods in addition to being optimal 
(or quadrature optimal) methods in the Hilbert space set up. 

The numerical results of the previous chapters show 
amply clearly the enhanced accuracy advantage of H-optiraal 
methods over the usual ones. 

The only disadvantage of H-optimal methods seems 
to be the storage requirement for the variable coefficients 
of these methods. However# for the present day computers# 
this much amoxant of storage is not really a problem. 
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APPENDIX 


DrSCUSBION AND SUGGESTIONS FOR FUTURE RESEARCH 

Ab we have ment i oned in the beginning, 'the present 
thesis may be regarded as a First e;:ploratory study of Sard's 
approach to optimal multistep methods' . 1 he findings by t»ow 

may be summarised as followss 

(a) by definition, the methods have been derived for 

? 2 ■■ 

certain spaces of functions (namely, H“ and L (C^)) as a 

whole. This entails an <at pr-i€>r-i assumption that we are dealing 
with equations whose solutions are from these spaces, (b) The 
formulae are characterized by their being inter polatory for 
certain lo(.al functions and their behavior (c) as r oo and 

(d) as step size h -♦ O have been determined . It is shown that 
(p) under the hypothesis of Picard's theorem the methods are 
convergent. The important point about the convergence here is 
that there is no restriction on the solutions (that it be a 
member of the space of the derivation of the method), (f) The 
coefficients of the methods, which are variable, however, are 
sl.own to differ from those of the parent usual methods by 
0(h), h -♦ 0. (g) Q(h'^''^S bounds ot> local truncation errors are 

established, for p--th order parent metliods in both implicit 
and explicit cases. Also, (h) bounds for discretization & (i) 
round off errors have been obtained which as far as thie order 
if, concerned para] lei those of tlie parent usual methods. (j) 



The numerical experiments on two dozen problems, conducted on 
a fixed step size basis, has given encourag ing results. 

The present work is limited by obvious considerations, 
not the least amongst them being due to the mist that 
surrounds the beginnings of any new approach. For instance the 
analytical derivation, coupled with the 0<h) closeness of the 
coefficients, strongly suggests an asymptotic expansion in 
powers of h which will not only result in clearing up 
questions related with various notions of stability and 
sharper asymptotic behavior of the numerical solutions, but 
will also facilitate a computation of the coefficients 
specially required in an adaptive implementation of the 
methods with stepsize control. 

Several aspects of the problem which we could not study 
concern: (a) various notions of stability, (b) stiff and 
singular equations, <c) system of equations, (d> initial and 
boundary value problems for higher order- equations, (e) 
special classes o.f differential equations, <f) any serious 
study of formulae i nt ei'po 1 atory for any specific functions 
that have been studied earlier in other context e.g. by 
Gautschi CAS], Lambert CAIZ], Lanczos CA14], Cash CA5], Corrol 
CA5], Iserles [All], and Neta and Ford CA15] etc, (Our 
motivation regarding this point came from the work of Chawla 
and Kaui C16-18], Kaul C>0] and Finney and Price, Jr. [24] as 
mentioned in the first chapter.), (g) spaces other than H^<C^) 



and L^(6 ) and <h) experiments with different values of r for 
r 

problems for which solutions are expected to be with a graded 
smoothness . 

We acknowledge that in order to take extra advantage of 
the accuracy that these optimal formulae deliver, it is 
critical that r be chosen with care , From an intuitive 
standpoint, as a general guideline. and L^<C^> with a 

smaller r would be more appropriate when the actual solution 
IS expected to be less smooth and a larger value of r for 
solutions with a higher smoothness. 

This is also in agreement with the following table in 
which we give the ratio, for the Problems ? and 4, of relative 
errors of optimal as compared with the usual methods in the 
case of explicit methods : 

Table_Al 

Ratio of end point <X|^ = 1.2) relative errors for problems ? & 4 


Table 

No. 


Multistep Method (EXPLICIT) 


M 

M^ 

M^ 

M^ 

M^ 

m'^ 

2 . 9 

10.373 

10.730 

11.803 

11.566 

17.049 

11 . 803 

3.9 

8.119 

7.344 

17.106 

8.119 

3.693 

10.882 

4.9 

9 . 970 

9.708 

12 . 490 

10.230 

64.728 

11.177 

5 . 9 

6 .997 

9.987 

9 .696 

6 .833 

2 . 989 

12.787 


We note that the singularity in Problem 5 is nearer the 

circle C than it is for Problem 4. and the methods perform 
r 



nearer . 


relatively much better when the singularity is 
Moreover, in the quadrature optimal cases the gains are 
I'elatively higher as compared to the optimal casses. 

The corresponding table in the implicit case is as 
foil ows : 

Table A2 


Ratio of end point (X|^=1.2) relative errors for problems 3 & 4 


Tabl e 

No. 


Multistep Method (IMPLICIT) 



M 

M^ 

M^ 

M^ 

m"^ M^ 



2.10 

1.957 

1.959 

1 .796 

2.145 

1.431 

2.724 

1 .760 

3.10 

0.823 

-0.386 

1.139 

-O'. 421 

0.913 ■ 

-0.144 

8.716 

4.10 

5.632 

4.039 

1 .780 

7.263 

1.162 

4.292 

2.270 

3.10 

0.143 

-0.279 

0.660 

0 .230 

0.744 

0.312 

-7.265 


In the above table, we observe that in the implicit case 
tfie performance of quadrature optimal methods (corresponding 
to Tables 2.10 and 4.10) in all cases is as stipulated. In the 
optimal cases (corresponding to Tables 3.10 and 5.10), 
however, there is no consistent pattern. This may be because 
for both the Problems 3 and 4 the errors are already quite 
small and are of both positive and negative signs. 

In this connection, a theoretical determination of r, 
unfortunately involves a (hitherto unavailable) knowledge 
about the degree of approximation 



E <y.X)= infC || y-g !| : g € X^, l! g I1 < X }. 

1 s r 


where <C^ ) or is the norm in X^ and 

II. IL is an appropriate norm for the solution space (e.g., I lyl L 

S ^ 

+ A II y' 11^. II . 11^ denoting the sup norm). This is because a 

tighter bound for the local truncation error is given by 

I T <y) I < inf C E (y;X) + X I1 T |1 }. 

n r nr 

In view of the above, in a practical situation a 

determination of r, right now, lies in the domain of 

ex per i men tat ion. 

In the thesis we have mostly concentrated on theoretical 


aspects and there too on only those for which we had something 
definite to say. Thus, there remains a lot of work to be done 
both in tfiooretical as well as practical directions. 
Concerning stability, for instance, we have been involved only 
with the weaker version when the characteristic polynomial of 
a’s is stable < i. e, its roots are on or inside the unit 
circle and those of modulus 1 are simple). Even though for the 
mere convergence it is enough (as follows from Theorem 1 
Chapter 6>, for stronger versions and for avoiding other type 
of numerical instabilities one has to work with only those 
optimal methods which correspond to usual methods possessing 
the appropriate stability. Thus. for instance, 
quadrature optimal methods with s 0. paras. t.c roots could 
ruin the relative stability of the solution, and therefore 



should not be applied for problems where this is important. 

The note at which we ended Chapter 6 is based on our 
numerical experience for a fixed step size only. when the 
coefficients have already been determined . The practical 
situation, indeed, is far from it . For, fixed stepsize 
numerical results say little about the overall potential of 
the formulae. To be of any practical use these formulae must 
be competitive in a variable-step, variable-order environment. 
This requires an efficient separate routine for a 
determination of the variable coefficients of the optimal 
methods. This extra work involves setting up the linear system 
of normal equations and solving it. (If at a step one chooses 
m non-zero a's and n b's, it involves (m+l+n)n function 
(kernel and its derivative, etc.) evaluations and about n^/ 5 
arithmetic operations at that step.) The- use of one or more 
parallel processors for this purpose has also to be 
i nvest i gated , 

In the thesis we have considered only the scalar case. 
However, in practice one comes across the vector case of 
systems of equations. in which, all the more, different 
components might possess different degrees of smoothness. In 
principle, the ideas of the thesis naturally extend to this 
case in a similar fashion (in about the same way as the proof 
of Picard’s theorem generalizes from the scalar to the vector 
case). Here the inner product for the product space would be 



(y .Z) = J: (y. , 2 ) 

j 1 1 H. ■ 

Where is a space of choice for the i-th component. Thus 

H.-s could be H^(C^) or L^C^). with the same or different 
r’s, or they might be some other spaces altogether. It 
may be noted here that if the a's and H. 's to each component 
are the same, the cofficients of the optimal formula coincide 
with those of the scalar case. However if the b’s are kept 
independent for each component the normal equations decouple 
with each other with respect to H^'s. The earlier discussion 
regarding the choice of r applies mxitatis mutandis to this 
case, making it clear that this product also needs a lot of 


future investigation, 

' A step size control for a predictor - corrector pair is 
made possible in our case by Lemma 6, Chapter 6. If the 
predictor - corrector usual formulae are respectively of p-th 
and (p+l)-th order, so are the related optimal formulae, by 
the lemma. Thus .following Cash CA3] and Shampine and Gorden 


CA21], we may use the standard strategy : 

With y<x^^j> as the exact solution, ^^+2 predicted 

value and y^^2 corrected value at ><^^^2 ' 

v<Xn.l> = 

and '''Xn+i> = V^+i + 

If 0<hP'^^) terms are neglected. TE = yn+2-yn+l 
estimate for the local truncation error for y^^j and it is 
this error which we will control. The scheme is based on local 



extrapolation, that is our error estimate is for the predicted 
solution ’ while it is the asymptotically more accurate 

solution Vp+j which is actually carried forward. This 
procedure, where the step control algorithm is based on an 
estimate of the local error in the predicted, rather than the 
corrected solution, is a common one. With TE as before in the 
scalar case (and putting TE = 11 ^p + 2~^n + l ^ vector 

case) and denoting the local accuracy required by s, the way 
in which the step size is controlled, when using p-th order 
formula, is as follows ; 

<i) If TE > e, halve h and recompute the solution at • 

<ii) I f «: > TE > , accept the solution and keep h fixed; 


< 1 1 i ) If 

TE 

< e/Z^^^ , accept 

the 

SO 1 ut 1 on 

and 

doub 1 e 

h. 

It 

may 

be remarked here 

that 

we do 

not 

know 

presently 

whether 

our optimal formulae 

possess 

the 

usua 1 

type of 

asymptot i c 

error estimates. 

and 

consequently a 

step size 


control with respect to the corrector is faced with 
difficulties. 

Finally a list of some additional and recent references 
on multistep methods, which could not find place in the main 
body of the thesis, follows: 
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1984. 
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